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1  Variational  convergence  of  bifunctions 


The  analysis  of  bivariate  functions,  or  more  briefly  bifunctions,  and  that  of  their  maxinf -points,  has 
played  in  key  role  in  “Variational  Analysis”  at  the  conceptual,  theoretical  and  computationally  levels. 
Simply  think  about  the  Lagrangians,  and  augmented  Lagrangians,  and  the  role  they  play  in  the  de¬ 
velopment  of  duality  theory  as  well  as  in  the  design  of  algorithmic  procedures  such  as  interior  point 
methods  and  Sequential  Quadratic  Programming  to  name  just  a  couple  of  examples.  The  same  applies 
to  Hamiltonians  associated  with  Calculus  of  Variations  and  Optimal  Control  problems  as  well  as  to  the 
study  of  minimax  functions  initiated  by  Danskin  [12]  and  further  elaborated  and  exploited  computa¬ 
tionally  by  Polak  [30],  Konnov  [27]  and  their  associates.  There  is  also  an  extensive  class  of  problems 
that  aren’t  traditionally  formulated  in  this  mode,  ‘find  a  maxinf-point  of  a  bifunction, ’  but  can  be  recast 
advantageously,  theoretically  and  computationally,  in  this  framework.  To  bring  to  the  fore  this  unifying 
relationship,  more  specifically  as  far  as  approximation  issues  are  concerned,  is  the  preeminent  aim  of 
this  article.  In  the  process,  we  also  touch  on  existence;  the  interdependence  between  the  tools  used  to 
obtain  existence  for  a  number  of  these  problems  was  already  brought  to  light  by  Aubin  and  Ekeland 
[7,  Chapter  6].  In  infinite  dimensional  spaces  maxinf-points  problems  are  related  to  the  Mountain  Pass 
Theorem,  minimal  surfaces  problems,  the  Ky-Fan  inequality  and  the  Shauder  fixed  point  theorem.  In 
this  article,  we  deal  with  the  finite  dimensional  case 

More  explicitly,  given  a  bifunction  F  :  C  x  D  — >  M,  we  are  interested  in  finding  a  point,  say  x  £  C, 
that  maximizes  with  respect  to  the  (first)  x- variable  the  infimum  of  F,  inf y^D  F(-,y),  with  respect  to 
the  (second)  y-variable.  We  refer  to  such  a  point  x  as  a  maxinf-point  and  write 

x  £  argmaxinf  A  =  argmaxinf  F  =  argrnax  (  inf  F(x,  y)). 

( x,y)£CxD  x£C  y&D 

Our  concerns  extend  to  approximate  maxinf-points 

for  e  >  0,  x£  is  an  approximate  maxinf-point  of  F  if  F(x£,  ■)  >  supinf  F  —  e, 

that  both  from  an  application’s  viewpoint  and  operationally  might,  conceivably,  be  more  important 
than  maxinf-points  themselves. 

In  some  particular  instances,  for  example  when  the  bifunction  is  concave-convex,  a  maxinf-point 
can  be  a  saddle  point  but,  in  general,  it’s  just  a  maxinf-point,  or  when  minimization  is  with  respect  to 
the  first  variable  and  the  maximization  of  F  is  with  respect  to  the  second  variable,  a  minsup-point.  To 
study  the  stability,  and  the  existence,  of  such  points,  and  the  sensitivity  of  their  associated  values,  one 
is  led  to  introduce  and  analyze  convergence  notion(s)  for  bifunctions  that  in  turn  will  guarantee  the 
convergence  of  either  their  saddle  points  [6,  4,  9],  or  just,  their  maxinf-points. 

This  paper  is  devoted  to  describe  and  analyze  this  collection  of  problems  that  includes  (geometric) 
variational  inequalities,  standard  and  set-valued  fixed  points  problems,  linear  and  nonlinear  comple¬ 
mentarity  problems,  Nash  equilibrium  for  non-cooperative  games  and  Walras  economic  equilibrium 
problems.  Each  one  of  these  examples  deserves,  eventually,  a  more  comprehensive  analysis  than  what 
is  appropriate  to  include  in  this  motivating  article;  we  essentially  limit  ourselves  to  illustrating  how  the 
general  approach  is  applicable  in  all  these  various  instances. 

Our  major  tool  is  the  notion  of  lopsided  convergence,  introduced  in  [5]  but  revisited,  and  more 
extensively  analyzed  in  [25],  so  that  a  wider  class  of  applications  could  be  handled.  In  [25],  the  accent 
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was  placed  on  finite-valued  bifunctions  defined  on  a  product  set  which  would  allow  us  to  clarify  the 
connection  with  the  the  earlier  work  [5]  for  the  extended  real- valued  functions.  The  still  slightly  more 
general  definition  of  lopsided  convergence  on  which  we  rely  in  this  article  completely  cuts  its  embryonic 
cord  with  the  extended-real  valued  framework  and  allows  us  to  analyze  approximation  results  for  an 
even  wider  class  of  variational  problems. 

In  §2,  we  characterize  variational  inequalities,  xed  points  problems,  complementarity  problems, 
Nash  equilibrium  and  Walras  economic  equilibrium  as  maxinf-problems.  In  §3,  we  refine  and  extend 
a  collection  of  results  about  lopsided  convergence  that  provide  the  main  instruments  to  analyze  the 
stability  properties  for  each  one  of  these  problems  as  detailed  in  §4-10. 

2  Examples 

In  order  to  set  up  the  pattern  of  our  approach,  let’s  construct  for  a  few  problems  in  this  class  the 
associated  bifunctions  and  show  that  their  solutions  are  exactly  the  maxinf-points  of  these  bifunctions. 

2.1  Linear  complementarity  problems 

The  linear  complementarity  problem  (LCP)  can  be  formulated  as  follows, 

find  z  such  that  0  <  z  T  Mz  +  q  =  w  >  0, 

where  Manx  n-matrix,  q  £  Mn  and  a  1  b  means  that  two  vectors  are  orthogonal.  Let’s  associate  the 
bifunction 

K(z ,  v)  =  (Mz  +  q,v  —  z)  defined  on  R”  x  R" 
with  this  complementarity  problem. 

2.1  Proposition  (maxinf-points  and  solutions  of  LCP).  z  solves  the  linear  complementarity  problem 
if  and  only  if 

z  £  argmax2>0  [inf„>o  K(z,  v )  ]  and  K(z,  •)  >  0. 

Proof.  Suppose  z  solve  the  LCP.  Since  inf„>o  K(x,  v )  =  — oo  unless  Mz  +  q  >  0,  a  condition  satisfied 
by  z,  v  =  0  yields  the  minimum.  Hence,  0  =  max2>o  —  (Mz  +  q ,  z)  and  this  maximum  is  attained  by  z 
with  K(z,v )  =  0  +  (Mz  +  q,v)  >  0  for  all  v  >  0.  On  the  other  hand,  when  z  is  a  maxinf-point  of  K 
on  R”  x  R™  and  K(z,  •)  >  0,  it  implies  (i)  z  >  0,  (ii)  K(z,  0)  >  0  which  yields  (Mz  +  q,z)  <  0,  and 
(iii)  0  <  —(Mz  +  q,  z)  +  inf„>o(M£  +  q,  v)  implies  Mz  +  q  >  0.  Combining  (ii)  and  (iii),  one  obtains 
(Mz  +  q,  z)  =  0,  and  thus  z  also  solves  the  linear  complementarity  problem.  □ 

Stability  analysis  of  the  solutions  to  linear  complementarity  problems  can  thus  be  undertaken  in 
terms  of  the  stability  properties  of  the  maxinf-points  of  the  corresponding  bifunctions.  More  specifically, 
given  a  sequence  of  linear  complementarity  problems, 

find  z  such  that  0  <  z  ±  Muz  +  qv  =  w  >  0 

where  Mu  and  qu  converge  ‘appropriately’  to  M  and  q,  one  can  define  a  sequence  of  approximating 
bifunctions 

Kv(z ,  v )  =  (Mvz  +  qu,v-z)  defined  on  ]Rn+  x  Mn+ 
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and  recast  the  convergence  of  the  solutions  of  the  approximating  linear  complementarity  problems  in 
terms  of  the  convergence  of  the  maxinf-points  of  the  bifunctions  Ku . 

A  similar  argument,  would  work  when  dealing  with  a  nonlinear  complementarity  problem  (NCP), 

find  z  such  that  0  <  z  _L  H(z)  >  0 

where  H  :  Mn  x  Mn  is  a  vector-valued  function  and  the  associated  bifunction  K  :  1R7]  x  ]R!+  is  de¬ 
fined  by  K(z,v )  =  ( H(z),v  —  z).  However,  one  can  also  view  such  a  problem,  as  well  as  the  linear 
complementarity  problem,  as  a  special  case  of  the  next  class  of  problems. 

2.2  Variational  inequalities  and  fixed  points 

Let’s  consider  the  following  variational  inequality  (V.I.):  find  u  £  C,  a  non-empty,  convex  subset  of 
Mn .  such  that 

u  £  C,  — G(u )  £  Nc(u) 

with 


•  G  a  function  from  C  into  Mn,  usually,  but  not  necessarily,  continuous, 

•  Nc(u )  =  {z  £  Mn  \(z,u-u)  <  0,  V-u  G  C)  the  normal  cone  to  C  at  u. 

With  {Cu  C  Mn,  v  £  JV}  a  sequence  of  convex  sets  converging  to  C  and  {Gu  :  Cu  — >•  Mn ,  v  £  LV} 
a  sequence  of  continuous  functions  converging  (appropriately)  to  G,  we  like  to  find  conditions  under 
which  one  can  refer  to  the  variational  inequalities:  find 

u  £  Cv  such  that  —  Gu  (it)  £  Nc»  (u) , 

as  approximating  V.I.’s.  In  particular,  one  would  like  to  be  able  to  assert  that  the  solutions  of  the 
approximating  variational  inequalities  converge  to  the  solution(s)  of  the  limit  one. 

As,  in  the  case  of  the  linear  complementarity  problem,  the  approach  that  we  follow  is  to  reformulate 
the  problem  in  the  following  terms:  Define  the  bifunctions 

K(u,v)  =  (G(u),v  —  u)  on  dom  K  =  CxC 


and,  for  v  £  LV, 

Ku(u,v)  =  (Gu(u),v  —  u)  on  dom Ku  =  Cv  x  Cv . 

and,  as  the  next  proposition  shows,  the  solutions  of  the  variational  inequality  can  be  identified  with  the 
maxinf-points  of  the  corresponding  bifunction. 

2.2  Proposition  (identifying  V. I. -solutions  as  maxinf-points).  Consider  the  V.I., 

u  £  C  so  that  — G(u )  £  Nc(u). 
and  the  associated  bifunction  K  :  C  x  C  M  with 

K(u,  v )  =  ( G(u),v  —  u). 
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Then,  u  is  a  solution  of  the  variational  inequality  if  and  only  if  it’s  a  maxinf-point  of  K,  i.e.,  u  G 
argmaxc  g  where  g{u)  =  inf„gc  [K(u,  v )  ]  on  C2,  and  K(u ,  •)  >  0. 

Proof.  Observe  that  u  €  C  is  a  solution  of  the  variational  inequality  then  —G(u)  €  Nc(u)  implies 
K(u,-)  >  0  that,  in  turn,  implies  g(u)  >  0.  On  the  other  hand,  by  definition,  for  all  u  G  C,  g(u)  < 
( G(u),u  —  u)  =  0.  Consequently,  g{u)  =  0  and  hence,  u  maximizes  g  on  C,  or  still,  it’s  a  maxinf-point 
of  Ah 

Conversely,  if  u  maximizes  g  on  C  and  K(u,  •)  >  0,  then  (—G(u),v  —  u)  <  0  for  all  v  G  C,  i.e.,  it’s 
a  solution  of  the  variational  inequality.  □ 

Of  course,  the  same  argument  applies  to  the  approximating  variational  inequalities  and  their  cor¬ 
responding  bifunctions.  Convergence  of  the  solution(s)  of  the  variational  inequalities  can  thus  be  for¬ 
mulated  in  terms  of  the  convergence  of  the  maxinf-points  of  the  bifunctions  Kv  to  the  maxinf-points  of 
K.  The  main  issue  will  be  to  identify  the  appropriate  convergence  notion  for  bifunctions  that  will  yield 
the  convergence  of  these  maxinf-points. 

To  see  how  the  nonlinear  complementarity  problem, 

find  z  such  that  0  <  z  T  H(z)  >  0, 

is  equivalent  to  the  variational  inequality 

find  z  >  0  so  that  (—H(z),v  —  z)  <  0,  Vu  >  0, 

can  be  argued  as  follows: 

•  if  z  solves  NCP  then  ( H(z),z )  =  0,  z  >  0  and  H(z)>  0  imply  ( H(z),z )  =  0  >  (—H(z),v)  for  all 
v  >  0,  i.e.,  z  solves  the  V.I., 

•  if  z  solves  V.I.,  z  >  0  and  one  has  (H(z),z)  >  (H(z),v)  for  all  v  >  0.  Hence  H(z)  >  0  since 
otherwise  there  is  no  solution  (— >  —  oo).  Also,  (H(z),0)  =  0  >  (H(z),z)  and  hence  H(z)  _L  z,  i.e., 
z  solves  the  NCP. 

Although  the  details  will  be  provided  in  §6  and  §9,  let’s  already  suggest  how  fixed  point  problems 
are  related  to  solving  variational  inequalities  and,  consequently,  fit  into  this  class  of  maxinf-problems. 
Given  G  a  continuous  function  from  a  convex  set  C  into  C  C  Mn,  finding  x  such  that  G{x)  =  x  can 
be  reformulated  as  solving  the  variational  inequality  (x  —  G(x))  G  Nc(x)  since,  in  particular,  this  last 
inclusion  requires  that  a  solution  must  satisfy  \x  —  G{x)  |2  =  0. 

2.3  Non-cooperative  Games 

Consider  a  game  with  a  finite  collection  of  players  A:  For  each  player  a  G  A,  Ca  C  Mn  denotes  the  set 
of  available  strategies.  The  choice  of  xa  G  Ca  yields  a  return  ra(xa,  x_a)  where  x_a  is  the  vector  of  the 
strategies  selected  by  the  remaining  players:  A  \  {a}; 

V  (x  G  A  :  Ua  :  C &  x  Ca  — ^  fFl,. 

a£A\a 

2The  function  g  has  been  extensively  studied  in  the  variational  inequalities  literature  under  the  name  of  gap  function, 
a  specific  instance  of  a  merit  function  for  a  variational  inequality,  cf.  [17,  §1.5.3] 
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The  strategies  x*  =  (x*,  a£i)  determine  a  Nash  equilibrium  point  of  this  game,  when 

Va  £  A  :  x*a  £  argmaxXagCa  ra(xa,  x!a); 

for  further  reference,  let’s  denote  this  game  by  =  {(Ca,  ra)  |  a  £  „4}. 

For  all  a  £  A,  {C'f,  v  €  iV}  a  sequence  of  sets  converging  to  Ca  and  {rf,^  £  iV}  a  sequence  of 
payoff  functions  converging  in  an  appropriate  sense  to  ra,  we  are  looking  for  conditions  that  will  allow 
us  to  assert  that  the  games  Qv  =  {(C^,rf)  |  a  £  *4}  approximate  Q,  i.e. ,  that  the  Nash  equilibrium 
points  of  the  games  Qv  approximate  those  of  the  game  Q. 

Our  approach  to  obtain  existence  and  continuity  results  again  relies  on  setting  up  a  bifunction  N, 
known  as  the  Nikaido-Isoda  bifunction,  defined  as  follows:  N  :  C  x  C  — >•  1R  where  C  =  n«e.4  an<^ 

N(x ,y)  =  ^  (^k^-a)  -ra{ya,x-a)y 

4. 

2.3  Proposition  (Nash  equilibrium  points  as  maxinf-points).  The  strategies  x*  =  (x*,  a  £  A),  with 
x*  £  Ca  for  all  a  £  A,  determine  a  Nash  equilibrium  point  if  and  only  if 

(r“(Xa’  X-a )  _  ra(ya,  x*-a) )  >  0  for  all  ya  £  Ca, 

a£A 

or  equivalently,  if  and  only  if 


x*  is  a  maxinf-point  of  N  and  inf^gc  N(x*,y)>0, 
where  x  =  (xa,  a  £  A),y  =  (ya,  a  £  .4.). 

Proof.  If  x*  =  (x*,  a  £  A)  is  a  Nash  equilibrium  point,  for  all  a  £  A, 

ra(x*a,x*_a) -ra(ya,x*_a)  >  0,  V ya  £  Ca, 

and  consequently  the  sum  over  a  £  A  must  also  be  nonnegative.  On  the  other  hand,  if 


(ra«^-a)  -ra(ya,x*_a)^j  >0  Vya  £  Ca, 


aeA 


it  implies,  in  particular,  that  given  any  player  a  £  A,  for  all  ya  £  Ca, 


ra{x*a,x*_a) -ra(ya,x*_a)+  (ua'(x*a,,x*_a,)  -  ua'(x*a,,x*_a,)^J  >  0. 


a'GA\a 

The  second  term  in  this  sum  is  0,  and  thus  x*  £  argrriax^,  ra(xa,x*_a ),  i.e.,  x*  =  (x*,  a  £  A)  is  a 
Nash  equilibrium  point. 

Turning  to  the  second  identity  involving  N,  observe  that  for  all  x  £  C,  inf,y  N{x,  y)  <  0.  Indeed,  if 
for  all  a  £  A,  xa  £  Ca, 


inf  N(x,y ) 
yeC 


Y]  ra(xa,  x_a)  -  sup  ra(ya,x_a) 
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Clearly, 


Va  €  A  :  ra(xa,x_a)  -  sup^  ra(ya,X-a )  <  0, 

and  hence  inf yN(x,y)  <  0.  But  since  for  x*  =  (x*,  a  £  A),  a  Nash  equilibrium  point,  inf yN(x*,y)  = 
N(x*,x*)  =  0,  it  follows  that 

x*  £  argmax^  [  inf  N(x,y)]. 
y 

Conversely,  if  x*  =  (x*,  a  £  A)  £  C  is  a  maxinf-point  of  N  such  that  inf yec  N(x*,y)  >  0,  it  means 

(ra(XI>  x*-a)  -  ra(yai  Z-j)  >  0  for  all  ya  £  Ca, 

a€A 

i.e. ,  x*  is  a  Nash  equilibrium  point  as  follows  from  our  first  assertion.  □ 

One  defines  similarly  the  Nikaido-Isoda  functions  Nu  associated  with  the  games  Qu ,  and  in  view  of 
Proposition  2.3,  xv  =  (x„,  a  £  A)  is  a  Nash  equilibrium  point  for  Qv  if  and  only  if  it’s  a  maxinf-point 
of  Nv. 

Here  again,  the  question  of  the  convergence  of  Nash  equilibrium  points  can  be  formulated  in  terms 
of  the  convergence  of  the  maxinf-points  of  the  (bivariate)  Nikaido-Isoda  functions  Nv  to  those  of  the 
function  N. 

2.4  Walras  economic  equilibrium 

Our  last  example  is  a  classical  equilibrium  problem  in  Economics.  Here,  we  deal  only  with  the  Pure 
Exchange  model  or  equivalently  a  Walras  barter  problem,  but  one  can  easily  extend  the  result  to  the 
case  when  the  economy  also  includes  producers  [3].  The  economy  is  described  by 

£  =  {( ua,Ca,ea ),  a€  A}, 

where 

A:  the  finite  set  of  agents; 

ea  £  lRn:  agent’s  a  £  A  endowment ,  a  bundle  of  goods  to  be  traded; 

Ca  C  JRV+  non-empty,  convex  set  identifying  agent’s  a  survival  set, 
ua  ■  Ca  — >  dR,  agent’s  a  utility  function. 

Trading  takes  place  at  a  per-unit  market  price  pj  for  good  j,  j  =  1, . . . ,  n.  The  bundle  of  goods  agent 
a  could  acquire  is  limited  by  the  budgetary  constraint  (p,x)  <  (p,ea).  It’s  assumed  that  agents  act  as 
utility  maximizers.  Thus,  given  p  £  JRf ,  each  agent  a  £  A  will  end  up  with  its  (consumption)  demand 

ca{p)  =  argrriaxxGC.;  { ua(x )  |  (p,x)  <  (p,ea)} 

which  we  assume  to  be  well  (uniquely)  defined.  Note  that  ca(p )  =  ca(ap)  for  any  positive  scalar  a,  i.e., 
the  agents  demand  functions  are  homogeneous  of  degree  0  with  respect  to  prices.  So,  we  may  as  well 
restrict  the  choice  of  p  to  A  =  {p  £  WCf  |  Pj  =  1 } ,  the  price  simplex. 

This  economy  £  is  operational  only  if  for  each  good,  total  supply  exceeds  total  demand,  i.e.,  if 

^  sa(p)  =  s(p)  >  0  where  sa(p)  =  ea  -  ca(j>); 

Ci^iA 
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this  is  the  market  clearing  condition.  The  function  s  :  — >•  Mn  is  called  the  excess  supply  function. 

Following  Walras  [36],  a  price  vector  p  £  A  so  that  s(p)  >  0  is  called  an  equilibrium  price  ( system ,)■,  the 
existence  of  such  an  equilibrium  price  isn’t  clear-cut,  albeit  well-known  since  the  mid-50’s  [13]. 

Insatiability ,  usually  expressed  in  terms  of  one  of  more  goods  being  persistently  attractive  to  some 
agents,  will  mean  that  at  equilibrium  whenever  si(p)  >  0,  i.e.,  good  l  is  in  excess  supply,  there  is  no 
interest  in  further  trading  and,  consequently,  pi  =  0.  On  the  other  hand,  when  pi  >  0,  indicating 
that  some  agents  might  still  be  interested  in  further  acquisition  of  this  good,  however  supply  has  now 
been  exhausted,  i.e.,  si(p)  =  0.  Thus,  every  equilibrium  price  system  satisfies  {s(p),p)  =  0  3.  Our 
predominant  objective  is  the  analysis  of  the  stability  of  equilibrium  prices  to  (global)  perturbations  of 
the  economy: 

-  perturbations  of  the  agents’  utilities, 

-  perturbations  of  the  agents’  endowments. 

For  this  purpose,  one  introduces  the  bifunction  W  :  A  x  A  — »  IR,  to  which  one  refers  as  the  Walrasian, 

W(p,q)  =  ( q,s(p )); 

and  the  ‘approximating’  Walrasians 

Wu(p.q )  =  (q,  sv(p))  on  A  x  A, 
where  su(p)  =  (e«  “  ca(P ))  and 

Ca(p)  =  argmaxx6Ca  «(x)  |  (p,x)  <  (p,eua)}. 

Eventually,  existence  and  stability  results  will  be  deduced  from  general  results  applied  to  the  functions 
W  and  Wu.  At  this  point,  our  only  aim  is  to  identify  the  maxinf-points  of  the  Walrasian  W  with  the 
equilibrium  prices  of  £. 

2.4  Proposition  (Walras  equilibrium  prices  and  maxinf-points).  Every  maxinf-point  p  of  the  Wal¬ 
rasian  such  that  W(p,  •)  >  0  on  A  is  an  equilibrium  point.  Moreover,  under  insatiability,  every  equilib¬ 
rium  point  p  of  £  is  a  maxinf-point  of  the  Walrasian  such  that  W(p,  •)  >  0  on  A. 

Proof.  Suppose  p  is  an  equilibrium  price  then  W(p,-)  >  0  on  A.  For  every  p  G  A  that’s  not  an 
equilibrium  point,  for  some  l,  si(p)  <  0,  i.e.,  inf q£A.{q,  s(p))  <  0  which  in  turns  implies  that  p  £ 
argmaxinf  W  since  for  equilibrium  point,  under  stationarity,  0  =  rnaxinf  W  =  {p,s(p)}.. 

Conversely,  if  p  is  a  maxinf-point  of  the  Walrasian  with  W(p,-)  >  0,  it  follows  that  for  all  unit 
vectors  eJ  =  (0, . . . ,  1, . . . ,  0),  the  jth  entry  is  1,  (e3 ,  s(p))  >  0  which  implies  s(p)  >0.  □ 


3  Lopsided  convergence 

Having  suggested  a  unified  framework  for  this  family  of  variational  problems,  one  can  reasonably  expect 
that  also  that  a  significant  number  of  their  properties  can  be  brought  to  light  when  examining  those  of 

3Jesse  Holzer  pointed  out  the  need  to  include  explicitly  the  insatiability  condition 
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the  maxinf-points  of  the  associated  bifunctions4.  That’s  certainly  the  case  as  far  as  existence  of  solutions 
is  concerned  and,  as  will  now  be  laid  out,  the  stability  of  these  solutions  under  various  perturbations. 
The  major  tool  in  deriving  these  stability  results  is  the  notion  of  lop-convergence,  or  lopsided  convergence 
for  bifunctions. 

In  this  section,  not  only  do  we  review  the  basic  definitions  but,  in  the  process,  expand  the  family 
of  bifunctions  to  which  the  basic  results  derived  in  [25]  are  applicable.  Precisely,  to  deal  with  these 
applications,  we  had  to  revisit  the  ‘classical’  definition  of  convergence  for  bifunctions.  The  (almost) 
classical  framework,  going  back  to  Rockafellar’s  work  on  saddle  functions  [32,  Sections  33-39],  was 
primed  in  considering  bifunctions  that  are  extended  real- valued  and  defined  everywhere  on  the  product 
of  two  linear  spaces,  here  Mn  and  Mm.  It  turns  out  that  this  paradigm  comes  with  unnecessary 
technical  obstacles  on  the  theoretical  side  and  renders  it  ungainly  when  dealing  with  the  rich  family 
of  applications  being  considered  here.  This  has  led  us  to  a  ‘parallel’  framework  where  the  bifunctions 
— and  also  functions  in  the  univariate  case —  are  real- valued  and  only  defined  on  the  product  of  certain 
subsets  of  Mn  and  Mm.  In  [25]  we  restricted  our  attention  to  a  class  of  bifunctions  that,  via  a  ‘natural’ 
embedding,  allowed  us  to  connect  lopsided  convergence  of  these  finite-valued  bifunctions  to  that  of 
‘proper’  extended-real  valued  bifunctions.  This  relationship  is  now  lost  and  it  doesn’t  seem  that  it  can 
be  recovered.  The  definition  of  lop-convergence  is  necessarily  one-sided,  in  that  one  is  either  interested 
in  the  convergence  of  maxinf-points  or  minsup-points  but  not  both;  in  general,  the  maxinf-points  are 
not  minsup-points  and  vice-versa.  Here,  as  in  [25],  definitions  and  results  are  stated  for  the  ‘maxinf’- 
case.  We  conclude  the  section  with  an  application  of  lopsided  convergence  to  obtain  an  extension  of  Ky 
Fan’s  Inequality  [19]  to  situations  when  the  domain  of  definition  is  not  necessarily  compact;  for  further 
analysis  refer  to  [25]. 

Henceforth,  the  term  bifunction  is  reserved  for  finite-valued  bivariate  functions  defined  on  the  prod¬ 
uct  of  two  non-empty  subsets  of  Mri  and  ttlm .  One  writes, 

fv-biv(Mn+m)  =  {F:  C  x  D  -A  I?  |  0  +  C  C  Mn ,  0  /  C  Rm} 
to  denote  this  class  of  bifunctions.  For  F  G  fv-biv(lRn+rn),  let 


supinf  F  denote  the  ‘optimal’  value,  i.e. ,  resulting  from  maximizing  on  C,  inf^g i>F(-,y), 
x  is  a  maxinf-point  of  F  G  fv-biv(Mn+rn)  if 

x  G  argmaxinfF  =  argmaxinf  F  =  argrnax  (  inf  F(x,y)j. 

(. x,y)eCxD  xGC  y&D 


and 


for  e  >  0,  xe  is  an  approximate  maxinf-point  of  F  if  F(xe ,  •)  >  supinf  F  —  e 


the  set  of  all  such  approximating  maxinf-points5  is  denoted  by  &argmaxinf  F. 

4  One  might  also  expect  that  the  design  of  algorithmic  procedures  can  also  be  guided  by  this  viewpoint  as  has  been  the 
case  in  our  own  work  on  stochastic  equilibrium  problems  [15,  14] 

Alternative  definitions  of  an  approximating  maxinf-point  could  be  considered,  for  example,  one  would  allow  also  for 
approximating  the  infimum  with  respect  to  y  as  well  as  the  supremum  with  respect  to  x.  Although,  this  gets  notionally  a 
bit  more  involved,  it  pretty  much  leads  to  same  type  of  results  than  those  developed  here.  In  fact,  in  many  applications, 
the  infimum  in  the  second  variable  is  actually  attained  and  a  maxinf-point  becomes  a  maxmin-point  however,  because  the 
bifunction  is  usually  not  concave  with  respect  to  the  first  variable  getting  to  such  a  maxmin-point,  often,  turns  out  to  be 
quite  challenging. 


3.1  Definition  (lop- convergence,  fv- biv).  A  sequence  in  fv-b\v{]Rn+m),  :  Cv  x  Dv  -> 
lop-converges,  or  converges  lopsided,  to  a  function  F  :  C  x  D  — »•  ]R,  also  in  fv-biv(lRn+Tn),  if 

(a)  for  all  y  £  D  and  all  ( xv  £  Cu)  -4i£C,  there  exists  €  Du)  — >  y  such  that 

limsup vFv{xv,yv)  <  F(x,y), 

(b)  for  all  x  £  C,  there  exists  (xv  £  Cv)  — >  x  such  that  given  any  (yu  £  Du)  — >  y, 

liminf„  Fv(xv,  yu)  >  F{x,  y)  when  y£D  and  Fv(xv,  yu)  — >  oo  when  y£D. 

Lop-convergence  is  ancillary  tight  when  (h)  is  strengthened  to 

(b-t)  (b)  holds  and  for  any  e  >  0  one  can  find  a  compact  set  Be,  possibly  depending  on  {xv  — >•  x}, 
such  that  for  all  v  sufficiently  large, 

infcm Be  Fv(xu ,  •)  <  info*  Fu(xu,  •)  +  e. 

Finally,  it’s  said  to  be  tight  if  it’s  ancillary-tight  and  (a)  is  strengthened  to 

(a-t)  (a)  holds  and  for  all  e  >  0  there  is  a  compact  set  As  such  that  for  all  v  large  enough, 

suPxeC-'nAe  inf y&D-  Fu(x,y)  >  supx6C,  miy&D,  Fu(x,y)  -  e, 

Not  only  does  the  following  results  apply  to  a  wider  family  of  bifunctions  but  they  sharpen  markedly 
the  assertions  of  [25,  Theorem  4]  As  will  be  made  clear  in  the  subsequent  sections,  they  also  allow  for 
substantial  improvements  of  the  quite  limited  collection  of  results,  available  in  the  literature,  when  they 
are  applied  to  our  family  of  variational  problems. 

3.2  Theorem  (convergence  of  maxinf-points) .  When  the  bifunctions  {Fu }^gjv  lop-converge  ancillary 
tightly  to  F,  all  in  fv-biv(fRn+m)  with  supinf  F  finite,  and  £u\£  >  0,  then  every  cluster  point  x  £  C 
of  a  sequence  of  £v-maxinf-points  of  the  bifunctions  Fu  is  a  £-maxinf-point  of  the  limit  function  F. 

In  particular,  this  implies  that  in  these  circumstances,  every  cluster  point  of  a  sequence  of  maxinf- 
points  of  the  bifunctions  Fu  is  a  maxinf-point  of  the  lop-limit  function  F. 

Proof.  Let  xv  £  £v-  argmaxinf  Fu  and,  without  loss  of  generality,  assume  that  the  entire  sequence 
xv  — >  x  £  C  and  let  ye  be  such  that  F(x,ye)  <  inf jjF(x,-)  +  £.  That  x  £  £- argmaxinf  F  is  inferred 
from  the  following  string  of  inequalities, 

F(x,  y£)  >  limsup,,  Fv[xv ,  yv)  >  liminf„  Fv [xv ,  yu)  >  liminf^ (supinf  Fu  —  £v)  >  supinf  F  —  e. 

The  first  one  of  the  inequalities  follows  from  3.2(a)  which  provides  an  appropriate  sequence  (yu  £  Du)  — > 
y£.  The  next  one  is  immediate.  The  third  one  follows  from  Fu(xu ,  yv)  >  inf pv  Fu(xu ,  •)  >  supinf  Fv — £v. 
The  last  one  relies  on  3.2(b-t)  as  is  shown  next. 

Indeed,  with  x  £  C  such  that  inf^  F(x,  •)  >  supinf  F  —  £,  let  (xu  £  Cu )  — >  x  be  one  of  the  sequences, 
and  Be,  the  corresponding  compact  set,  postulated  by  3.2(b-t),  such  that  for  v  sufficiently  large 

(i)  V  (yv  £  Du )  — >  y,  liminf uFu(xu,yu)  >  F(x,y)  or  Fu(xu ,yu)  F oo  depending  on  y  £  D  or  not, 

(ii)  By"  £  Du  n  Be  such  that  Fu(xL',yu)  <  inf^i-  Fv(xu ,  ■)  +  £u]  recall  £u\£. 
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Passing  to  a  subsequence,  if  necessary,  let  yu  — >  y.  Either  Fu(xu,  yv)  /'oo  in  which  case  so  does 
inf £)v  Fu(xu,  •)  and  also  supinf  Fu  which  certainly  implies  lim  inf^(supinf  Fu  —  £u)  =  oo  >  supinf  F  —  e. 
Or,  liminf„  Fu(xu,  yv)  <  oo,  and  one  has 

liminfj, (supinf  Fu  —  eu)  >  liminfj,  ( \nijjv  Fv{xv ,  ■)  —  e„)  >  ( lim  inf  Fu(xLI ,  yv)  —  ej) 

>  F(x,  y)  —  e  >  inf£>  F{xu ,  •)  —  e  >  supinf  F  —  e 

again  yielding  lim  inf „  supinf  [Fv  —  ev)  >  supinf  F  —  e.  □ 

When  lopsided  convergence  is  (fully)  tight,  some  additional  properties  turn  out  to  be  valuable. 

3.3  Theorem  (under  tight  lop-convergence).  When  the  bifunctions  {Fv}v^]\f  lop-converge  tightly  and 
supinf  F  is  finite,  then  supinf  Fu  — >  supinf  F .  In  fact,  if  for  any  e„  \0  and  xv  £  e^-argmaxinf  Fu ,  then 
Ejrsupinf  Fu  =  infer'  Fv(xv ,  ■)  — >  supinf  F . 

Proof.  It  suffices  to  show  that  limsup^  supinf  Fu  <  supinf  F  since  we  just  finished  showing  that  under 
ancillary  tightness,  lim  inf  v  supinf  Fu  >  supinf  F]  setting  £v  =  e  =  0.  When  the  sequence  {Fu,  v  G  IV} 
is  (fully)  tight,  in  view  of  3.2(a-t),  given  s  >  0,  for  all  u  large  enough,  one  can  always  find  xv  £  (Cu rAe) 
such  that  inf^i-  Fu(xv ,  •)  >  supinf  Fv — e  with  Ae  the  pertinent  compact  set.  Passing  to  a  subsequence,  if 
needed,  xv  — >  x  in  Ae.  Let  y  £  D  be  such  that  F(x,  y)  <  inf/j  F(x,  ■)+£  <  supinf  +£  and  (yu  £  Du)  — >•  y 
such  that  lirnsup u  Fv (xu ,yu)  <  F(x,y).  One  has 

£  +  supinf  F  >  F(x,  y)  >  lirnsup, y  Fv(xl/ ,  yu)  >  limsup^(inf£)^  Fu(xu ,  •)  >  limsup„  supinf  Fu . 

Since  this  holds  for  all  e  >  0,  supinf  Fu  — >  supinf  F  again  follows. 

The  same  argument  applies  when  supinf  Fv  is  replaced  by  e^supinf  Fu  with  £v  \  0  and  rather  than 
calculating  sup^c^  inf'/>  Fu(x,  •)  one  is  content  with  a  e^-maxinf  point.  □ 


3.4  Corollary  (existence  of  approximating  solutions).  When  the  bifunctions  {Fu}u g^v  lop-converge 
tightly  and  supinf  F  is  Unite.  If,  for  epsilon  >  0,  x  is  a  £-maxinf-point  of  F,  one  can  always  find 
sequences  £u\£  and  xp  £  e^-maxinf  Fu  such  that  xv  — >  x.  Conversely,  if  for  e  =  0  such  sequences  exist, 
then  supx  miy  Fu{x,  y)  — >  infy  F(x,  y). 

Proof.  We  already  proved  that  under  these  hypotheses  supinf  Fv  — >  supinf  F  which  also  implies  that 
for  v  sufficiently  large,  supinf  Fu  is  finite-valued.  Thus,  the  functions  x  i— >  inf^  Fv[x,  ■)  are  well- 
defined  finite  valued  on  a  nonempty  domain.  We  obtain  the  hypo-convergence  of  these  functions  to 
x  i— >  'mipF^x,  ■)  from  [25,  Theorem  3].  The  assertions  then  follow  from  [25,  Theorem  2].  □ 

It  might  be  useful  at  this  point  to  emphasize  the  role  played  by  ancillary  tightness,  cf.  Theorem  3.2.  Let’s 
begin  with  a  simple  example  where  Fv  lop-converge  to  F ,  supinf  F  is  finite  and  argmaxinf  F  ^  0  but 
supinf  Fv  74  supinf  F  precisely  because  for  all  £  >  0  there  is  no  compact  set  Be  — and  a  corresponding 
index  u£ —  such  that  for  x  £  Cu ,  inf b£  Fu(x,  •)  <  inf  Fu(x,  ■)  +  £  for  all  v  >  ve.  Let  Cv  =  C  =  1R  and 


Fu{x,y) 


0  when  y  ^  {0,  u} 
<  1  for  y  =  0 
^  v  for  y  =  v. 
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Clearly,  the  Fv  lop-converge  to  F  for  F(x,  0)  =  1  and  F(x,y)  =  0  when  y  /  0.  supinf  F  =  1  is  finite 
and  attained  at  any  point  in  M  x  {0},  but  supinf  Fu F oo  ^  l!  It’s  obvious  that  there  is  no  compact 
set  and  an  index  set  with  the  desired  properties.  This  example  would  also  work  defining  Fu(x,  v)  =  k, 
with  k,  any  finite  number  greater  than  1  instead  of  Fu(x ,  v)  =  v. 

In  particular,  this  implies  that  the  existence  of  an  maxinf-points  doesn’t  allows  us  to  conclude  that 
the  ‘required’  compact  sets  exist.  One  can  translate  this  as  follows:  the  existence  of  a  maxinf-point  for 
the  limit  problem  doesn’t  guarantee  ancillary  tightness,  as  one  might  have  expected  or  hoped  for.  The 
example  shows  that  one  must  be  concerned  about  other  maxinf-points  that  disappear  at  (  ~  converge 
to)  the  horizon. 


3.5  Definition  (Ky  Fan  functions).  A  bifunction  F:  C  x  C  — >•  1R  in  fv-biv(M2n)  with  C  a  non-empty 
convex  subset  of  ]Rn  such  that 

(a)  V  y  G  C:  x  >->•  F(x,  y)  is  use  (upper  semicontinuous)  on  C, 

(b)  Vr  €  C:  y  F(x,y)  is  convex  on  C. 
is  said  to  be  a  Ky  Fan  function. 

Note  that  the  set  C  is  not  required  to  be  compact.  However, 

3.6  Lemma  (Ky  Fan’s  Inequality;  [19],  [7,  Theorem  6.3.5]).  Suppose  F  :  C  x  C  — >•  1R  is  a  Ky  Fan 
function  with  C  compact  and  such  that  F(x,x)  >  0  (on  C  x  C).  Then,  the  set  of  maxinf-points  of  F 
is  a  nonempty  subset  of  C .  Moreover,  for  every  maxinf-point  x  of  C,  F(x ,  ■)  >  0  on  C . 

When,  the  domain  is  not  compact,  one  can  nevertheless  obtain  the  existence  of  maxinf-points  by  relying 
on  [25,  Theorem  8]  which  shows  that  the  lop-limit  of  a  sequence  of  Ky  Fan  functions  is  also  a  Ky  Fan 
function.  Consequently, 

3.7  Corollary  (extension  of  Ky  Fan’s  Inequality;  [25,  Theorem  9]).  If  F  is  a  Ky  Fan  function  defined 
on  C  x  C  with  0  ^  C  convex  and  one  can  find  sequences  of  compact  convex  sets  j Cu  C  Mn}  and 
(finite-valued)  Ky  Fan  functions  {Fv :  Cu  x  Cv  — >■  lR}u<ziv  lop-converging  ancillary  tightly  to  F,  then 
every  cluster  point  x  of  any  sequence  {xv ,  v  £  IN}  of  maxinf-points  of  the  Fu  is  a  maxinf-point  of  F . 

3.8  Remark  (a  word  of  caution).  It  should  be  kept  in  mind  that  all  one  can  expect  from  this  general 
approach  when  applied  to  specific  instances  is  that,  mostly,  one  can  come  up  with  sufficient  conditions. 
For  example,  in  §6  when  dealing  with  the  convergence  of  fixed  points,  we  place  ourselves  in  a  ‘standard’ 
environment,  i.e.,  the  limit  problem  is  to  find  x  in  a  compact,  convex  set  C  C  lRn  that  is  a  fixed  point 
of  a  continuous  function  G  :  C  — »  C  and  this  limit  problem  is  being  approximated  by  a  sequence  of 
problems  having  similar  characteristics.  Of  course,  one  can  trivially  find  examples  of  mappings  G ,  not 
necessarily  continuous,  that  map  C  into  C  but  nonetheless  have  a  fixed  point  and  the  approximating 
sequence  would  just  consist  of  problems  reproducing  the  limit  one.  In  such  a  situation,  one  might  even 
be  able  to  prove  lop-convergence,  but  the  conditions  of  Theorem  6.2  wouldn’t  apply.  The  same  type  of 
observation  could  be  made  in  the  case  of  variational  inequalities,  cf.  §5,  as  well  as  to  any  other  class  of 
variational  problems  analyzed  here. 
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4  Complementarity  Problems 

Let’s  now  return  to  the  linear  complementarity  problem, 

find  z  >  0  so  that  Mz  +  q  >  0,  (Mz  +  q)  _L  z, 

and  a  sequence  of  ‘approximating’  (truncated)  linear  complementarity  problems, 

find  z  €  [0,  rv]  so  that  Mv z  +  qu  >  0,  (Muz  +  qv)  _L  z, 

where  Mv  — >•  M,  qu  — »•  q  and  for  j  =  1, . . . ,  n,  0  <  rj  — >■  oo  as  v  — >  oo.  We  associate  the  bifunction 

IF(z,u)  =  (Mz  +  q,v-z),  K:  x  ^  ->  1?, 

with  the  linear  complementarity  problem  and  the  bifunctions 

Kv(z,  v )  =  (ATz  +  g",  v  -  z),  FT :  [0,  r"]  x  IR\  -)•  1? 

with  the  ‘approximating’  problems.  We  used  the  word  ‘approximating’  with  quotes,  because  these  are 
genuine  approximating  problems  only  if  it  can  be  shown  that  these  approximating  problems  actually 
generate  approximating  solutions  of  the  given  linear  complementarity  problem.  In  view  of  our  analysis 
in  §2.1,  in  particular  Proposition  2.1,  and  the  results  in  §3,  one  could  validate  the  term  approximating 
by  finding  conditions  under  which  the  bifunctions  Ku  lop-converge  to  K. 

4. 1  Theorem  (lopsided  convergence  of  LCP) .  As  long  as  Mv  — >•  M,  qu  — >•  q  and  rv  /■  oo,  the  bifunctions 
Ku  lop-converge  to  K.  If  in  addition,  P  =  {z  >  0  |  Mz  +  q  >  0}  is  included  in  the  inner  limit, 
Liminf^oo  Pu,  of  the  polyhedral  sets  Pu  =  {z  G  [0  ,ru]  \  Muz+qu  >  0},  then  the  sequence  {Ku ,  v  G  IN} 
lop-converges  ancillary  tightly  to  K  which  means  that  any  cluster  point  of  the  solutions  of  the  sequence 
of  ‘truncated’  linear  complementarity  problems  is  a  solution  of  the  limiting  linear  complementarity 
problems. 

In  fact,  this  is  also  the  case  if  rather  than  maxinf-points  of  the  approximating  problems,  one  would 
allow  for  just  s^-maxinf-points  with  £j,\0. 

Proof.  For  any  sequence  {zu  £  [0,  ru]}uejy  converging  to  z  >  0  and  any  v  >  0,  ( Muz u  +  qu,  v  —  zv)  — > 
(Mz  +  q,v  —  z)  yields  condition  3.1(a);  the  case  zv  — >  z  with  x  ^  iR”  doesn’t  have  to  be  considered 
since  no  sequence  {zv  €  [0, ry]}^ejv  can  converge  to  a  point  z  ^  J?".  Finally,  given  z  €  If",  eventually 
z  G  [0,r"]  for  all  v  sufficiently  large,  say  v  >  vz ,  since  the  entries  of  the  vector  rv  all  converge  to  oo. 
So,  given  z  >  0,  define  the  sequence  zv  — >  z  to  be  any  point  zv  G  [0,  rv]  for  u  <  vz  and  zv  =  z  when 
v  >  vz.  Clearly,  such  a  sequence  converges  to  z  and  whatever  be  the  sequence  vv  >  0  — >  v  >  0,  one  has 
(Mv zv  +  qv ,vv)  — >  (Mz  +  q,v),  i.e.,  condition  3.1(b)  is  also  satisfied. 

Ancillary  tightness  requires  in  addition  that  given  any  sequence  zv  G  [0,  rv\  z,  for  all  e  >  0,  one 
can  hnd  a  compact  set  Be  and  an  index  ve  such  that  for  v  >  v£, 

inft,GRjnss  (Mvzv  +  qu ,v)  <  inf v<zRn(Mv zv  +  qv,v)  +  e; 

note  that  Be  isn’t  necessarily  a  rectangle  of  the  type  [0,  ru\ .  This  condition  can  only  be  satisfied  if  the 
right-hand  sides  stay  bounded,  this  requires  that  the  (limiting)  linear  program, 

min(Mz  +  q,  v)  so  that  v  >  0, 
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and  the  approximating  ones,  for  v  sufficiently  large, 

min (Muzv  +  qu,v)  so  that  v  >  0, 

be  bounded,  i.e. ,  their  optimal  solutions  is  v  =  0  with  optimal  value  0.  This  will  only  occur  if  Mz  +  q 
and  Mv z  +  qv ,  for  v  sufficiently  large,  are  both  non-negative.  So,  to  satisfy  condition  3.1(b-t),  it’s 
necessary  that  for  all  z  £  P  =  { z  G  M+  \  Mz  +  q  >  0 } ,  one  can  find  a  sequence  zv  £  Pu  =  { z  £ 
[0,^]  |  Mvz  +  qv  >  0}  converging  to  z.  This  means  precisely  that  P  must  be  included  in  the  inner  limit 
of  the  polyhedral  sets  Pu ,  cf.  [33,  Definition  4.1].  □ 

One  can  seek  conditions  that  will  guarantee:  the  inner-limit  of  the  polyhedral  sets  Pv  is  included 
in  P.  As  it  turns  out,  as  long  as  the  sets  Pu  are  non-empty  for  v  sufficiently  large,  P  always  includes 
the  outer  limit,  lim  sup^^^  Pv ,  when  Mu  — >  M,  qu  — >•  q  and  the  entries  of  ru  — »•  oo.  This  follows 
almost  immediately  from  the  definition  [37,  Proposition  1],  Thus,  we  are  essentially  requiring  that 
P  =  lim^oo  Pu ,  i.e.,  P  is  actually  the  limit  of  the  polyhedral  sets  Pv .  A  substantial  literature, 
surveyed  and  complemented  in  [37,  28],  has  been  devoted  to  this  issue.  Hence,  we  won’t  deal  here  with 
all  specific  instances  that  are  of  particular  interest  in  the  theory  underlying  the  linear  complementarity 
problem,  cf.  [11],  that  would  lead  us  too  far  astray  from  the  main  theme  of  this  paper.  Let’s  just  record 
a  couple  of  sufficient  conditions  that  will  serve  as  examples. 

Let’s  suppose  the  approximating  problems  are  formulated  without  upper  bounds  on  z,  i.e., 

find  z  €  J?,”  so  that  Muz  +  qu  >  0,  (Muz  +  qu )  T  z, 
with  associated  bifunctions 


Kv(z,v)  =  ( Muz  +  qu,v-z),  Ku :  1R”  x  lRn+ ^  JR. 


Let  posA  denote  the  positive  hull  of  the  columns  of  the  matrix  A.  Then,  ancillary  tightness  of  the 
collection  of  bifunctions  [Ku ,  v  £  LV}  can  be  characterized  as  follows: 


4.2  Corollary  (continuity  of  polyhedral  set-valued  mappings).  With  P  =  {z  >  0  |  Mz  +  q  >  0}, 
Pu  =  {z  >  0  |  Muz  +  qv  >  0},  then  P  =  Lim^oo  Pv  if  and  only  if 


pos 


i  0  I  T  • 

ID  Limsup  pos 

U  /  v — ^OO 


(Mu)t 

~qv 


I  0 
0  -1 


where  q  and  qv  are  the  row-versions  of  these  vectors.  Thus,  under  this  last  condition,  the  bifunctions 
Ku  are  ancillary  tight,  and  consequently,  any  cluster  point  of  solutions  of  the  approximating  LCP’s  is 
a  solution  of  the  limiting  linear  complementarity  problem. 

Proof.  It  suffices  to  appeal  to  [37,  Theorem  3].6. 


References  [35,  28]  provide  a  number  of  sufficient  conditions.  Situations  that  are  more  amenable  to 
immediate  verification  can  be  found  in  [37],  for  example: 

6for  a  slightly  sharper  but  more  involved  condition,  cf.  [28,  Corollary  4.6] 
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4.3  Corollary  (non-empty  interior  criterion,  [37,  Corollary  7]).  Suppose  that  for  v  sufficiently  large, 
the  polyhedral  set  Pu  and  P  have  non-empty  interior  and  no  row  of  the  matrix  [M,  q ]  is  identically  0, 
then  Pu  P. 

Of  course,  one  could  consider  other  approximation  schemes  than  the  one  discusses  so  far.  One  could 
rely  on  truncations:  for  example,  for  a  sequence  0  <  ru  Poo, 

find  z  €  [0,  ru]  so  that  Mu z  +  qu  €  [0,  ru],  ( Muz  +  qu)  -1  z, 

with  corresponding  bifunctions 

K£(z,v)  =  ( Muz  +  qu,v-z ),  Kvn:  [0,rv]  x  [0 ,rv]->R. 

Since  [0.ru]  is  compact  convex,  it’s  known  that  problems  of  this  type  always  have  a  solution,  see  [24], 
for  example. 

4.4  Theorem  (lop-convergence  of  LCP,  variant).  The  bifunctions  lop-converge  to  K  when  Mv  — >• 
M,  qu  — >•  q  and  ru  Poo.  If  in  addition,  Liminf u  P"  D  P  =  JR”  O  {z  |  Mz  +  q  >  0}  where 

P"  =  [0,rv}p{Z\Mvz  +  qv  >  0}, 

then  lop- convergence  is  ancillary  tight,  and  this  means  that  any  cluster  point  of  a  sequence  of  solutions 
of  the  approximating  problems  is  a  solution  of  the  (given)  limit  problem. 

Proof.  The  same  arguments  as  those  used  to  obtain  lop-convergence  and  ancillary  tightness  in  Theorem 
4.1  also  work  here;  except  that  now  the  sequence  vv  — >  v  is  such  that  vu  €  [0,  ru].  □ 

One  could  also  rely  on  the  truncation  used  by  Gowda  and  Pang  [21],  and  more  recently  by  Flores-Bazan 
and  Lopez  [20],  in  their  stability  analysis  of  the  solutions  of  linear  complementarity  problems.  Namely, 
let  0  <  d  £  Mn  and  av  Poo,  a  sequence  of  positive  scalars.  The  approximating  problems  are 

find  2  6  =  {z  €  J?"  |  ( d ,  z)  <  au]  so  that  Mv z  +  qv  >  0,  (Muz  +  qv)  _L  z, 

with 

P£  =  A1'  n{z\M"z  +  qu  >  0,  } 

and  associated  bifunctions 

Kua(z,v )  =  {Muz  +  qu,v-z),  Kua:  AAA"4  JR. 

Again,  since  Au  is  compact  convex,  it’s  known  that  problems  of  this  type  always  have  at  least  one 
solution. 

4.5  Theorem  (lop-convergence  of  LCP,  another  variant).  The  bifunctions  KA  lop-converge  to  K  when 
Mu  — >  M,  qu  — >■  q  and  au  Poo.  If  in  addition,  Liminf VPA  D  P  =  Av  fi  {z  |  Mz  +  q  >  0}  then  lop- 
convergence  is  ancillary  tight,  and  it  means  that  any  cluster  point  of  a  sequence  of  solutions  of  the 
approximating  problems  is  a  solution  of  the  (given)  limit  problem. 
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Proof.  Similar  to  that  proof  of  Theorem  4.4  except  that  the  sequence  vu  — >•  v  now  has  vv  £  A v .  □ 

Existence  of  solutions  to  the  truncated  LCP  problems  is  well-known.  Here,  one  could  derive  it 
directly  from  Ky  Fan’s  Inequality  3.6  since  the  sets  [0,  ru],  as  well  as  Au ,  are  non-empty,  compact 
and  convex  and  the  functions  K □  and  Kv^  are  Ky  Fan  functions,  cf.  §3.  But  even  ancillary  tight  lop- 
convergence  of  these  functions  doesn’t  settle  the  question  of  the  existence  of  solutions  to  limit  LCP. 
This  can  be  illustrated  by  the  following  simple  example. 


4.6  Example  (ancillary  tightness  and  existence  of  solutions).  Let 


Mv 


'0  u~r 
1  1  ’ 


0 

1  1  ’ 


then  Pp  — »  P  =  0  and,  consequently,  the  lop-convergence  of  the  bifunctions  K ^  to  K  is  ancillary 
tight,  but  clearly  the  solutions  zv  =  (zz,  0)  of  the  truncated  problem  don’t  converge  to  a  solution  of  the 
(limit)  linear  complementarity  problem  that  happens  to  have  no  solution;  Gowda  and  Pang  [21]  and 
Flores-Bazan  and  Lopez  [20]  go  through  a  detailed  analysis  that  allows  then  to  guarantee  the  existence 
of  solutions  for  matrices  M  that  fall  in  very  specific  classes. 

Deriving  simple  verifiable  conditions  that  enable  us  to  assert  that  the  sequence  is  tight  requires  a 
challenging  study  of  quadratic  forms  that  goes  beyond  the  scope  of  this  article. 


5  Convergence  of  variational  inequalities 

We  return  to  variational  inequalities  as  introduced  in  §2.1,  i.e.,  let  0  ^  C  C  Mn ,  convex,  compact  and 
G  :  C  —>  fRn  a,  not  necessarily  continuous,  function.  The  problem  is 

find  u  £  C  such  that  —  G(u)  £  Nc(u) 

where  Nq(u)  is  the  normal  cone  of  C  at  u.  Such  a  problem  can  be  reformulated  (Proposition  2.2)  as 
finding  the  maxinf-point  of  the  following  Ky  Fan  bifunction,  namely, 

K(u,v )  :=  ( G(u),v  —  u)  with  domK  =  C  x  C. 

Clearly,  this  bifunction  is  convex  (linear)  in  v,  it’s  use  in  u  on  C  when  G  is  continuous  and  moreover, 
K(u,u )  =  0.  Ky  Fan’s  inequality  3.6  then  guarantees  the  existence  of  maxinf-point,  say  u,  that  is  a 
solution  of  this  variational  inequality  with  K(u ,  •)  >  0. 

The  variational  inequality  problem  can  also  be  written  as 

find  u  £  C  such  that  ( G(u),v  —  u)  >  0,  Vu  £  C. 

In  this  context,  given  e  >  0,  the  concept  of  a  e-approximate  solution  ue  means  that  the  inequalities 
{G(ue),v  —  ue)  >  0  might  fail,  for  some  v  £  C  but  never  by  more  than  e.  Implicitly,  the  problem 
becomes 

find  u  £  C  such  that  ( G(u),v  —  u)  >  —e,  Vu  £  C; 
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it’s  essentially  this  approach  that’s  used  in  the  design  of  a  solution  procedure  for  stochastic  variational 
inequalities  in  [10]. 

From  the  convergence  results  in  §3,  Corollary  3.7  and  Theorem  3.2,  it  follows  that  a  sequence  of  Ky 
Fan  bifunctions  {Ku}u€n  lop-converges  ancillary  tightly  to  the  bifunction  K  will  yield  a  maxinf-point 
of  the  limit  bifunction  K  if  the  sequence  of  maxinf-points  of  the  bifunctions  Kv  admits  a  cluster  point. 
The  implication  for  variational  inequalities  is  the  following: 

5.1  Proposition  (convergence  of  variational  inequalities).  The  variational  inequality:  find  u  £  C  such 
that  — G(u )  £  Nq(u )  where  G  :  C  —>■  JRn  and  0  /  C  C  fRn  is  convex  will  have  a  solution  as  long  as  we 
can  find  a  sequence  of  sets  {Cu}u£]y  convex,  compact,  converging  to  C,  and  a  sequence  of  lsc  functions 
{Gv :  Cv  -3  fRn}u(z]v  converging  continuously  to  G  with  respect  to  the  sequence  Cu  -3  C. 

Proof.  The  functions  Gu  :  JRn  -3  JET1  converge  continuously  to  G  (on  the  sequence  Cv  -3  C ),  i.e.,  for 
all  xv  £  Cv  -3  x  €  C,  Gv{xv )  -3  G{x).  Define 

Ku(u,x)  =  (Gv(u),v -u)  on  dom  Ku  =  CuxCu. 

Then,  Ku  are  Ky  Fan  bifunctions  that  converge  lopsided  ancillary  tightly  to  K,  i.e.,  any  cluster  point 
of  the  solutions  of  the  approximating  variational  inequalities  is  a  solution  of  the  limit  one.  □ 

Requiring  continuous  convergence  might  appear  to  be  bit  too  stringent,  but  possibly  only  in  some 
very  specific  instances,  it’s  unavoidable.  To  support  this  assertion,  let’s  go  through  two  completely 
different  arguments.  One  that  studies  the  convergence  of  variational  inequalities  in  terms  of  generalized 
equations,  or  equivalently,  as  set-valued  inclusions.  A  second  one  that  consider  the  special  case  when 
the  variational  inequalities  can  be  identified  with  (convex)  optimization  problems. 

5.2  Example  We  consider  the  same  collection  of  variational  inequalities  with  Cv  — >  C  and  Gv  con¬ 
verging  appropriately  to  the  continuous  function  G  adding,  already,  the  assumpion  that  C  is  compact. 
It’s  this  ‘ appropriately ’  that  we  want  to  investigate  rewriting  the  variational  inequalities  as  set-valued 
inclusions, 

u  £  C  such  that  A(u )  3  0  where  A(u)  =  Nc(u)  +  G(u) 

and  for  u  =  1, . . . , 


u  £  Cv  such  that  Au(u)  3  0  where  Au(u)  =  Nc^(u)  +  Gu(u) 

Detail.  The  assumption  that  Cu  -3  C  and  C  convex  and  bounded  means  that  for  v  sufficiently  large 
enough,  all  sets  Cv  lie  in  a  bounded  region  [33,  Corollary  4.12].  Without  loss  of  generality,  we  may  as  well 
assume  that  all  sets  Cu  C  U  a  bounded  subset  of  JRn '.  Taking  this  for  granted,  we  now  concern  ourselves 
with  conditions  under  which  the  solutions  of  Au(u)  3  0  will  converge  to  those  of  A(u)  3  0.  In  view  of  [33, 
Theorem  5.37],  the  question  essentially  boils  down  to  asking  if  Au  A.  Now,  Cu  -3  C,  by  Attouch’s 
Theorem  [33,  Theorem  12.35]  yields  Nc^  ^3  Nc,  in  fact  because  these  are  convex  cones  Nc»  -A-  Nq, 
i.e.  totally  converge  to  [33,  Theorem  4.25(b)].  By  continuity  rgeG,  the  range  of  G,  is  bounded  in  C 
and  gph  G  is  clearly  connected,  it  follows  that  for  v  sufficiently  large,  gph  Gu  are  uniformly  bounded  if 
Gu  ft '3  G,  cf.  again  [33,  Corollary  4.12],  Moreover,  in  this  situation,  when  gphG'1'  -3  gphG,  they  actually 
totally  converge  [33,  Theorem  4.25(d)].  There  now  remains  only  to  appeal  to  [33,  Exercise  4.29(c)]  to 
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conclude  that  Au  ^  A  when  Gu  and  Nc  D  (gph  —G)°°  =  {0}.  This  last  condition  is  innocuous  since 
G  is  bounded  on  C.  We  now  turn  to  [33,  Corollary  5.45]  to  conclude  that  since  G  is  single-valued  and 
bounded,  then  for  all  practical  purposes  Gv  G  is  equivalent  to  Gv  converges  continuously  to  G.  □ 

Let’s  note  that  Robinson  [31,  22],  working  to  obtain  error  bounds,  under  significantly  more  restrictive 
conditions  (subdifferentiability  type-conditions)  is  also  led  to  require  continuous  convergence  of  the 
mappings  Gv  to  G.  In  [23,  §3],  Giirkan  and  Pang  make  a  thorough  analysis  of  the  convergence  of 
the  variational  inequalities  associated  with  finding  Nash  equilibrium  points  of  non-cooperative  games. 
Their  results  are  a  bit  more  specific  since  they  deal  with  a  class  of  non-cooperative  games  where  the 
strategies  sets  are  ‘fixed’,  i.e. ,  don’t  change  with  the  convergence  parameter  u.  Interestingly  enough, 
they  are  also  led  to  impose  continuous  convergence  on  the  gradients  of  the  reward  functions  that  in  this 
framework  correspond  to  the  continuous  convergence  of  the  functions  the  Gu . 

One  might  still  hope  that  one  could  escape  ‘continuous  convergence’  in  the  absolute  nicest  of  all 
situations. 

5.3  Example  We  again  consider  the  same  collection  of  variational  inequalities  with  Cu  — »  C  and 
Gu  converging  appropriately  to  the  continuous  function  G  but  adding  now  the  assumption  that  the 
functions  Gu  and  G  are  monotone  with  domain  lRn  . 

Detail.  The  functions  Gv  can  be  viewed  as  the  gradients  of  convex  functions,  say  gu  and  g,  and  the 
solutions  of  the  variational  inequalities  are  the  optimal  solutions  of  the  convex  optimization  problems 

min  gu  +  ic"  and  min  g  +  tc 

where  ld  is  the  indicator  function  of  the  set  D.  The  question  about  the  convergence  of  the  solutions 
then  conies  down  to  the  epi-convergence  of  these  functions  gv  to  g.  Since  Cu  — >  C ,  ic "  ic,  the 
sums  will  converge  when  dorng  and  C  cannot  be  separated,  certainly  satisfied  when  doing  is  Mn ,  and 
the  functions  gv  epi-converge  to  g  [33,  Exercise  7.47(b)].  But,  now  again,  via  Attouch’s  Theorem  [33, 
Theorem  12.35],  gu  -£->■  g  cannot  occur  unless  their  gradients  Gu  =  'Vgl/  converge  graphically  to  G  =  V g 
which  brings  us  back  to  continuous  convergence  via  [33,  Theorem  12.35]  already  cited  earlier.  □ 

Let’s  conclude  this  section  by  summarizing  our  results  as  follows: 

5.4  Theorem  (convergence  of  the  solution  sets  of  V.I.).  Suppose  j Cv  C  lRn ,  u  G  LV}  is  a  sequence 
of  convex,  compact  sets  converging  to  C  ^  0,  necessarily  convex  but  not  not  necessarily  compact,  and 
the  functions  {Gu  :  lRn  -3  lRn,  v  G  iV}  converge  continuously  to  a  (continuous)  bounded  function  G. 
Then,  the  solution  set 

D  =  {u\G(u)  +  Nc(u)  3  0} 

of  the  limiting  variational  inequality  contains  the  outer-limit  of  the  solutions  sets  of  the  variational 
inequalities,  i.e.,  Limsup^  Du  C  D  where 

Du  =  {u  \  Gu(u)  +  NC"(u)  3  0},  v  £  JV. 

Moreover,  if  u  €  D,  then  there  exists  approximate  solutions  uv  G  Cu  of  the  variational  inequalities 
— Gv(u )  G  Nc^(u)  3  0  such  that  uu  -3  u. 
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Proof.  Let’s  begin  by  observing  that  in  these  circumstances  the  sets  Du  are  nonempty  as  observed  at 
the  beginning  of  this  section  and  (uu  G  Du )  — >  u  yields  u  G  D  by  Proposition  5.1. 

On  the  other  hand,  if  u  G  D,  then  u  G  argmax(inf„ec(G(ri),  v  —  u))  which  in  view  of  Corollary  3.4 

uec 

means  that  there  exists  {eu\0,uL'  €  ev-  argmaxn€Cl,  [m.iv&cv  {Gu  {u),v  —  u)]},  or  equivalently 
and  (uv  €  Cu)  — >  u  such  that  (—Gv(uv),v  —  uv)<£v,  V v  €  Cv , 
which  yields  the  sequence  of  approximating  solutions.  □ 


6  Convergence  of  fixed  points 

Brouwer’s  Fixed  Point  Theorem,  and  its  classical  generalizations,  can  be  derived  from  Ky  Fan’s  in¬ 
equality.  Indeed,  with  C  a  non-empty,  compact,  convex  subset  of  lRn  and  G  :  G  — >■  G  a  continuous 
function,  let’s  define  the  bifunction 

F(x,  y)  =  (x  —  G(x),y  —  x),  i.e. ,  F  :  C  x  C  — >  ]R. 

F  is  clearly  a  Ky  Fan  bifunction  defined  on  a  product  of  compact,  convex  sets  and  for  all  x  G  Mn , 
F(x,x)  =  (x  —  G(x),x  —  x)  >  0.  Hence,  from  Lemma  3.6  follows  the  existence  of  a  maxinf-point  x  of 
F  such  that 

Vy  G  IB  :  F(x ,  y)  =  (x  —  G(x),y  —  x)  >  0. 

Since  G(x)  G  G,  recall  G  :  G  — >•  G,  and  since  F(x,  ■)  >  0,  again  by  Lemma  3.6,  choosing  y  =  G(x),  one 
has 

(x  —  G(x),  G(x)  —  x)  =  —  \x  —  G(x)|2  >  0 

and  this  can  only  occur  when  G{x)  =  x.  We  encapsulate  this  conclusion  in  the  well-known  theorem' : 

6.1  Theorem  (existence  of  a  fixed  point).  Let  G  :  G  — >•  G  be  continuous  where  C  C  Mn  is  (nonempty), 
compact  and  convex.  Then,  G  admits  a  fixed  point,  i.e.,  for  some  x  G  G,  G(x)  =  x. 

Let’s  now  turn  to  approximation  issues  and  consider  the  following  situation,  the  compact,  convex  sets 
Cu  — >•  G  which  is  nonempty,  convex  and  compact.  In  particular,  this  means  [33,  Corollary  4.11]  that 
for  v  sufficiently  large,  all  sets  Cv  are  contained  in  G  +  ylB  for  some  rj  >  0.  The  continuous  functions 
Gu  :  Cu  — >•  Cu  are  converging  continuously  to  G  :  G  — >•  G  with  respect  to  Cv  — >•  G:  V  xv  G  Cu  — >  x  G  G, 
Gu(xu)  — >■  G(x).  It’s  then  straightforward  to  verify  that  the  bifunctions  Ku(x,y )  =  (x  —  Gv(x),y  —  x) 
defined  on  G^  x  G"  lop-converge  tightly  to  K  :  (x  —  G(x),  y  —  x)  defined  on  G  x  G  and,  via  Theorems 
3.3  and  3.2,  this  implies, 

6.2  Theorem  (convergence  of  fixed  points).  In  the  situation  described  here  above, 

1.  for  all  v,  the  functions  Gv  have  a  least  one  fixed  point  in  Cv  as  is  also  the  case  for  G  in  C , 

2.  for  C #  the  set  of  cluster  points  of  all  the  fixed  points  of  the  functions  Gu ,  0  /  G#  C  G, 

3.  if  x  G  G#,  then  x  is  a  fixed  point  of  G  on  C , 

7in  [5]  a  coercivity  condition  had  been  imposed  on  G  :  (x,  G( x)}  <  |x|2  to  obtain  this  result  via  the  Ky  Fan  inequality. 
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4.  if  x  G  C  is  a  fixed  point  of  G,  there  exist  approximate  fixed  points  of  Gu ,  xv  G  Cu  converging  to  x. 

5.  if  x  is  a  cluster  point  of  Ev-inaxinf-points  of  the  functions  Fv{x ,  y)  =  (x  —  G(x),y  —  x)  defined  on 
Cu  x  Cv  with  £u\0,  then  x  is  a  fixed  point  of  G  on  C. 

Proof.  All  the  assertions  are  covered  by  the  paragraph  preceding  the  statement  of  the  theorem.  □ 

7  Stability  and  existence  of  Nash  equilibrium 

We  next  turn  to  non-cooperative  games  as  introduced  in  §2.3  and  deal  with  the  existence  and  the 
stability  of  Nash  equilibrium  points.  We  consider  a  game  Q  =  {(Ca,  ra)  |  a  G  .4}  and  an  approximating 
sequence  Qv  =  {(C%,  rf)  |  a  €  A}  with  a  finite  number  \A\  of  players.  For  each  player  a  £  A,  the  sets  Ca 
or  C'f,  subsets  of  Bn,  determine  the  available  strategies  xa  and  the  associated  reward  is  ra(xa,  x_a),  or 
rf(xa,  x_a),  where  x_a  is  the  vector  of  the  strategies  selected  by  the  remaining  players.  Nash  equilibrium 
points  are  strategies  that  satisfy 

M a  G  A:  x*  G  argmaxXaeCa  ra(xa,  x*_a)  or  x G  argmax^g^  ra(xa,  xv_fa ). 

Existence  and  continuity  results  for  Nash  equilibrium  points  will  be  derived  via  the  properties  of  the 
maxinf-points  of  the  bivariate  Nikaido-Isoda  bifunctions,  N  :  C  x  C  — »  1R  where  C  =  n„«&, 

N(x,y)  =  Yi  hot xa,x_a)  —  ra(ya,x_a)j , 

cl£.vA. 

and  Nu  :CU  x  Cu  ->■  B  where  Cv  =  UaeA  C ", 

Nu(x,y)  =  (<(*o>*-o)  -ra(ya,x-a)y 

a£A 

e-approximate  solutions  are  of  a  totally  different  character  than  those  for  variational  inequalities  or 
fixed  point  problems.  Here  it  means  that  if  we  add  up  by  how  much  each  individual  player  fails  to  meet 
its  maximum  rewards  the  total  won’t  turn  out  to  be  no  more  than  e;  when  evenly  split,  each  individual 
player  might  fall  short  by  no  more  than  e/|*4.|. 

Existence,  a  well-known  result  (cf.  [7,  Theorem  4.2],  for  example),  is  obtained  here  as  a  direct 
consequence  of  Ky  Fan  Inequality’s  3.6,  and  is  extended  by  appealing  to  Corollary  3.7;  refer  to  the  last 
few  paragraphs  of  §3  and  the  fact  that  the  Nikaido-Isoda  bifunction  N  is  finite-valued  on  C  x  C. 

7.1  Theorem  (existence  of  Nash  equilibrium  points).  If  for  all  a  G  A  the  sets  Ca  are  convex  and 
compact  and  the  Nikaido-Isoda  bifunction  N  :  Bnx ^  x  Bnx ^  — >  B  satisfy: 

(a)  for  all  y  G  Bnx  x  >->■  N{x,y)  is  use, 

(b)  for  all  x  G  Bnx^,  y  *->•  N(x,y)  is  convex, 
then  Q  has  a  Nash  equilibrium  point. 

Proof.  One  simply  appeals  to  Ky  Fan’s  Inequality  3.6  after  observing  that  N(x,x )  >0.  □ 

It’s  immediate  from  the  definition  of  the  Nikaido-Isoda  bifunction  that  it  will  be  a  Ky  Fan  bifunction 
under  the  following  conditions: 
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7.2  Proposition  (Nikaido-Isoda  as  a  Ky  Fan  bifunction).  A  sufficient  condition  for  N  to  be  a  Ky  Fan 
bifunction  is  the  following:  for  all  a  €  A,  the  sets  Ca  are  convex  and, 

a)  ra  is  use  and  for  all  xa  £  IRn ,  ra(xa,  •)  is  lsc; 

b)  for  all  x  £  IRnx^,  ra(-,x_a)  is  concave. 

Let’s  now  turn  to  stability  issues  related  to  perturbations  of  both  the  strategy  sets  and  the  payoffs.  For 
this  purpose  we  introduce  the  following  convergence  notion  for  a  sequence  of  (approximating)  games: 
=  {(Cf,rf),a  £  A}  for  ueW. 

7.3  Definition  (convergence  of  non-cooperative  games).  Convergence  of  a  sequence  of  games  {Gu ,  v£ 
IV }  to  a  game  Q  is  defined  in  the  following  terms:  for  all  a  &  A, 

a)  the  nonempty  compact  convex  sets  C f  converge  to  the  nonempty  compact  set  Ca; 

b)  the  sequence  of  Nikaido-Isoda  bifunctions  Nu  associated  with  the  games  Qv  lop-converge  ancillary 
tightly  to  the  Nikaido-Isoda  bifunction  N  associated  with  game  Q . 

When  the  collection  of  sets  Cva,  nonempty,  compact,  convex  converge  to  Ca,  this  limit  set  is  also  con¬ 
vex  [33,  Proposition  4.15].  Also,  if  Cf  — >  Ca  for  all  a  £  A,  then  ^ a  ~ ^  riae.4  [33,  Exercise  4.29]. 

An  extension  of  continuous  convergence  usually  defined  relative  to  a  fixed  set,  say  C,  cf.  [33,  §7.C] 
turns  out  to  be  a  sufficient  condition  for  lopsided  convergence  of  the  Nikaido-Isoda  bifunction  sequence 
Nu .  Indeed,  a  sequence  of  functions  {fv  :  Mn  — >•  1R,  v  £  IN}  is  said  to  converge  continuously  to 
f  :  IRn  — >  1R  relative  to  a  sequence  of  sets  Cu  — >  C  if  for  all  xv  — >  x  such  that  for  all  v  £  IV,  xv  €  Cv: 
fu(xu)  —>■  f(x).  Thus,  it’s  easy  to  see  that  condition  a)  in  the  definition  of  convergence  of  games  jointly 
with  the  continuous  convergence  of  rf  relative  to  Cf  for  all  a  £  A  imply  the  ancillary  tight  lopsided 
convergence  of  Nikaido-Isoda  bifunctions  associated  with  the  games  Qv . 

Combining  the  previous  results  and  observations,  we  can  formulate  our  first  stability  result  for  Nash 
equilibrium  points. 

7.4  Theorem  (convergence  of  Nash  equilibrium  points).  Suppose  the  games  {Gu ,  v  £  IV  }  converge 
to  a  game  Q  and  the  Nikaido-Isoda  bifunctions  Nv  associated  with  them  are  Ky  Fan,  then,  there  exist 
strategies  {xu  =  (xf,  a  £  A),  v  £  IN}  that,  for  all  v  £  IN,  are  Nash  equilibrium  points  of  Gu.  Moreover, 
any  cluster  point,  say  {x  =  (xa,  a  £  .A)}  of  such  a  sequence  of  Nash  equilibrium  points  is  itself  a  Nash 
equilibrium  point  of  Q. 

Proof.  Since  the  games  \GV ,  v  £  IN }  converge  to  a  game  G  and  the  bifunctions  Nu  are  finite  on 
nae^  Cf,  we  obtain  from  Theorem  7.1  and  the  properties  of  lopsided  convergence,  cf.  Theorem  3.2,  the 
assertions  of  the  theorem.  □ 

The  next  corollary  translates  the  stability  result  in  terms  of  the  original  formulation  of  the  games. 

7.5  Corollary  (sufficient  conditions  for  convergence  of  Nash  points).  Suppose  the  games  \  GV ,  v  £  IN} 
are  such  that  for  all  a  £  A  the  sets  Ca  are  convex  and  compact  and  that  the  payoff  functions  satisfy: 

a)  rf  is  use  and  for  all  xa  £  IRn ,  rf(xa,  •)  is  lsc; 

b)  for  all  x  £  Mnx^,  rf(-,x_a)  is  concave. 

Suppose  also  that  for  all  a  £  A,  rf  continuously  converges  relatively  to  Cf.  Then,  there  exist  a  Nash 
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equilibrium  {xu  =  (x^,  a  £  A),  v  £  IV}  of  Qv  for  all  u  £  IV,  and  any  cluster  point  of  this  sequence  is  a 
Nash  equilibrium  of  the  game  Q. 

Independently,  and  at  about  the  same  time  we  started  to  circulate  some  of  the  results  in  this  article, 
Giirkan  and  Pang  [23]  developed  an  alternative  approach  to  the  convergence  of  Nash  equilibrium  points 
based  on  multi- epi  convergence ,  more  precisely,  multi-hypo  convergence  since  the  formulation  is  in  terms 
of  rewards-maximization  for  the  players.  The  sequence  of  games  Qv ,  they  consider,  is  less  general  than 
those  being  analyzed  here.  In  their  framework,  the  strategies  sets  Ca  are  constant,  i.e.,  don’t  depend 
on  the  approximation  parameter  v.  We  extend  Gtirkan  and  Pang  [23]  definition  and  allow  for  this 
dependence  in  the  definition  below.  Let 


C-a  =  IL  Co!  and  similarly  Cv_a 


n 


cv 

a'e.4\{a 


7.6  Definition  (multi-hypo  convergence  a  la  Gurkan-Pang) .  For  all  a  £  A,  suppose  that  for  all 
{xu_a  £  Cf_a,  v  £  IV }  converging  to  X-a,  the  functions  {r^(-,  xf_a),  v  £  IV }  hypo-converge  to  ra(-,X-a); 
for  the  definition  of  hypo-convergence  for  finite-valued  functions  defined  on  subsets  of  Mn ,  cf.  [25,  §2], 

We  now  relate  the  (extended)  Giirkan-Pang  approach  to  that  based  on  lopsided  convergence. 


7.7  Theorem  (hypo-convergence  of  the  reward  functions).  Suppose  the  games  {Gu,  v  £  IV  }  are  such 
that  their  reward  functions  multi-hypo  converge  to  the  reward  functions  of  Q  =  {(Ca,ra),a  £  A}. 
Then  the  corresponding  sequence  of  Nikaido-Isoda  bifunctions  {Nu,  v  £  IV}  lop-converges  to  N,  the 
Nikaido-Isoda  bifunction  associated  vith  the  game  Q. 

Proof.  For  each  player  a  £  A,  hypo-convergence  of  the  functions  {r£(-,  xf_a),  v  £  IV  }  to  ra(-,x-a) 
means  [25,  Proposition  2] 

a) -a°°)  Mxu  £  C"  — >  x  £  Ca,  limsup„  r^{xu)  >  ra(x)  and  r^(aT)\  —  oo  when  x  Ca, 

b)  \/xa  £  Ca,3xa  — >  xa  such  that  liminf ur"(xu)  >  ra(x). 

Let’s  first  verify  that  condition  3.1(a)  is  satisfied,  namely  for  all  y  £  C  and  xu  £  Cu  — >•  x,  3  ( yv  £ 
Cu)  — >  y  such  that  limsupj,  Nu(xu,  yv)  <  N(x,y)  when  x  £  C  and  limsup„  Nu(xu,  yu)  \  —  oo  when 
x  £  C,  or  equivalently  that  given  any  y  £  C  and  any  xu  £  Cu  — >•  x  £  C  one  can  find  yu  —>  y  such  that 

limsupj,  {ra{xu)  -  ra(yva,x-aav ))  <  Y  (r“(x)  “  ra(ya,X-a)) 


when  x  £  C  and  otherwise 


limsupj,  Y  (ra(xu)  -  ra{ya,x-aau))  \  -  oo. 

CL£l>A. 

Condition  b)  guarantees  that  for  all  ya  £  C  and  xv_a  £  Cfa  — >  X-a  £  C-a,  one  can  find  y"  £ 
Cq  — y  ya  £  Ca  such  that  liminf^  r^(y^,  xv_a)  >  ra(ya,X-a).  On  the  other  hand,  given  the  sequence 
xu_a  £  Cf_a  — >  X-a  £  C-a  and  any  sequence  xva  £  Cva  either  limsup^  ra{xv)  <  ra(x  =  (xa,X-a))  or 
ra{xu)  \  —  oo.  This  means  that  for  any  sequence  xu  €  Cu  — >•  x  and  y  £  C  with,  for  all  a  €  A,  the 
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appropriate  choice  of  y%  —>■  ya  one  has 

limsup[r£(®,/)  -  rva(yva,xv_a)]  <  lim  sup  r£ (sc")  -  lim inf  rva (yva,  xv_a) 

V  V  v 

j  ra(x)  -  ra(ya,  x-a)  if  xa  G  Ca, 

-  [-00  if  Xa  £  Ca. 

Clearly,  this  remains  valid  once  we  take  the  sum  with  respect  to  a  G  A  and  thus  the  first  condition  for 
the  lopsided  convergence  of  the  sequence  Nikaido-Isoda  Nv  to  N  is  satisfied. 

Let’s  now  turn  to  verifying  3.1(b).  Given  x  G  C,  let’s  choose  our  sequence  (xva,xv_a)  G  Cv  — > 
(. xa,x-a )  in  such  a  way  that  for  each  a  £  A,  lim  inf ur^(xu)  >  ra(x)  as  predicated  by  the  hypo- 
convergence  of  the  reward  functions,  more  specifically,  by  assumption  b).  Given  any  sequence  {yu  G 
Cu ,  v  G  IN}  converging  to  y,  from  a)-a°°)  one  knowns  that  limsup„  r"(y",  xv_a)  <  ra(ya,  X-a)  when 
ya  G  Ca  or  rva{yva,xv_a)\  -  oo  otherwise,  or  equivalently,  liming  -rva{yva,  xv_a)  >  -ra(ya,X-a)  when 
Ua  G  Ca  or  r"(y",  xv_a)  /'oo  otherwise.  This  immediately  carries  over  to  the  Nikaido-Isoda  bifunctions 
by  taking  sums.  □ 

In  [23],  Giirkan  and  Pang  make  a  thorough  analysis  of  the  convergence  of  the  variational  inequality 
associated  with  finding  Nash  equilibrium  points  for  non-cooperative  games.  Their  results  are  a  bit  more 
specific  since  they  deal  with  a  class  of  non-cooperative  games  where  the  strategy  sets  are  ‘fixed’,  i.e., 
don’t  change  with  the  convergence  parameter  v.  Interestingly,  they  are  also  led  to  impose  continuous 
convergence  on  the  gradients  of  the  reward  functions  that  in  this  case  correspond  to  the  functions  Gv  in 
our  formulation  of  approximating  variational  inequalities  in  §5;  refer  also  to  the  more  detailed  comments 
about  continuous  convergence  in  that  section. 

8  Existence  and  stability  of  Walras  equilibria 

One  can  rewrite  the  Walras  barter  problem  as  a  non-cooperative  game  and  then  apply  the  results  of 
the  two  preceding  sections,  in  particular  §7.  Let’s  begin  by  following  this  path  but  in  the  end,  as  we 
shall  see,  a  more  direct  approach  via  the  Walrasian  bifunction  turns  out  to  be  more  expedient. 

Returning  to  the  Walras  model  introduced  in  §2.4,  our  collection  of  players  will  consist  of  the  individual 
agents  a  G  A  in  addition  to  a  so-called  Walrasian  auctioneer,  a  panoptic  player  whose  main  function  is 
to  choose  a  market  price  system  aimed  at  securing  a  market  equilibrium.  The  reward  functions  are 

Va  G  A  :  ra(x,x-a,p)  =  ua(x)  when  (p,  ea  -  x)  >  0,  x  G  Ca 

and  for  the  Walrasian  auctioneer, 

rw(Pi  xa)  =  inf  (q,  Y]  _  (ea  -  xa))  if  0  <  p  ^  0 

q>  0  L — 'a£A 

where  xa  consists  of  the  strategies  of  all  the  individual  agents.  In  this  model, 

•  ra  doesn’t  depend  on  x_a,  it’s  included  in  the  arguments  simply  for  consistency, 

•  when  0  <  p  /  0,  f'wiP^A)  =  ~°°  unless  ^ag^(ea  —  xa)  >  0  in  which  case  rw  =  0,  its  upper 
bound. 
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Now,  suppose  that  (xA,p),  with  p  ^  0,  is  an  equilibrium  point  of  the  Walras  barter  model  which  implies 
that  the  excess  supply  J2a€A^ea  ~  —  0-  Then,  clearly,  for  all  a  £  A,  xa  £  argmax^,  ra(x,  X-aip)  and 

p  £  argmax0<p^0  rw(-,xX)-  The  reward  function  of  the  Walrasian  auctioneer  doesn’t  depend  explicitly 
on  p  but  actually  does  so  indirectly;  in  fact,  once  a  price  system  has  be  chosen  so  that  YlaeA(ea~ Xa)  — 
given  that  for  all  a  £  A,  xa  £  argmaxx  ra(x,  x_a,p),  in  theory  any  other  non-negative  price  system  would 
“maximize”  rw  but  that  wouldn’t  guarantee  that  each  xa  would  then  maximize  ra(-,x_a,p).  Hence, 
{xaiP)  is  an  Nash  equilibrium  point  for  the  non-cooperative  game  defined  by  the  rewards  functions 
((ra,a€  A),rw). 

On  the  other  hand,  if  (xa ,P)  is  a  Nash  equilibrium  point  of  a  game  with  reward  functions  ((ra,  a  £ 
A),  rw^ ,  then  necessarily  for  all  a  £  A,  xa  £  argma xx{ua(x)  |  (p,  ea  —  x  >  0}  and  J2aeA(ea  ~  —  0 

since  otherwise  argrnax  rw  would  be  empty.  This  means  that  0  <  p  ^  0  with  xA  is  an  equilibrium  point 
of  the  Walras  barter  model. 

One  can  then  proceed  to  writing  down  the  corresponding  Nikaido-Isoda  bifunction, 

N(((x,p),  (y,q))  =  ^2aeA  ( ra(xa,x_a,p )  -  ra(ya,  x_a,p))  +  (rw(p,xA)  -  rw(q,xA)) 

but  an  explicit  expression,  in  terms  of  the  given  utility  functions  and  the  associated  budgetary  con¬ 
straints,  gets  a  little  unwieldy  that  makes  the  analysis  and,  in  particular,  the  study  of  the  convergence 
properties,  more  involved  than  they  should  be. 

So,  let’s  proceed  as  in  §2.4  and  work  with  the  Walrasian  bifunctions  and  deduce  both  existence  and 
stability  from  the  general  results  for  bifunctions.  Recall  that  £  =  {(ua,  Ca,  ea),  a  £  A}  provides  the 
description  of  the  economy,  W(p,q )  =  (s(p),q)  on  A  x  A  is  the  Walrasian  with  A  the  unit  simplex  in 
Rn  and  s(-)  is  the  excess  supply  function. 

8.1  Theorem  Suppose  that  p  i— >  s(p)  is  use  on  A,  then  £  has  at  least  one  Walras  equilibrium  point, 
say  p.  Moreover,  W(p,-)  >  0. 

Proof.  When  s  is  use,  so  is  W(-,q)  for  all  q  £  A  and  hence  W  is  then  a  Ky  Fan  bifunction,  finite¬ 
valued  on  the  compact  convex  set  A  x  A.  Since  for  any  p  £  A,  W(p,p)  >  0,  Ky  Fan’s  inequality  3.6 
immediately  yields  the  existence  of  a  maxinf-point.  Moreover,  W(p,-)  >  0  when  p  is  a  maxinf-point, 
again  by  Lemma  3.6.  □ 

Conditions  on  the  original  data  of  the  economy  £  under  which  the  excess  supply  function  is  use  are 
provided  in  the  next  proposition;  basically,  the  same  conditions  as  those  used  by  Arrow  and  Debreu  [3] 
to  derive  their  existence  result. 

8.2  Proposition  If  for  all  a  £  A  the  utility  functions  ua  are  use  and  concave  on  Ca,  and  the  initial 
endowments  ea  £  int  Ca,  then  W(.,q)  is  use  on  A,  for  all  q,  i.e.,  W  is  a  Ky  Fan  bifunction. 

Proof.  This  is  an  immediate  consequence  of  the  classical  results  for  the  sup-projection  of  a  bivariate 
function,  here  sa(p)  =  ea  —  sup{u(x)  |  (p,  x  —  ea)  <  0,  x  £  Ca},  cf.  for  example,  [33,  Chapter  7].  □ 

Next,  we  consider  a  sequence  of  economies  £u  =  {(uua,  C%,  eva)  |  a  £  A },  that  can  be  interpreted  as 
perturbation  of  the  utility  functions,  survival  sets  and  initial  endowments  e^intC'^'. 
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8.3  Definition  A  sequence  of  economies  {Sv,  v  G  IV  }  is  said  to  be  converging  to  an  economy  8  if  the 
Walrasians  Wu(p,  q)  associated  with  the  economy  8V  lop-converge  ancillary  tightly  to  the  Walrasian  W 
associated  with  8. 

This  means  that  convergence  of  the  economies  is  defined  in  term  of  the  convergence  of  their  associated 
Walrasians.  Sufficient  conditions  to  guarantee  the  lop-convergence  ancillary  tightly  of  the  Walrasian 
bifunctions  Wv  were  given  in  one  of  our  earlier  paper  [26].  Indeed,  if  for  all  a  G  A  the  utility  functions 
v8a  are  use  and  concave  with  the  same  domain,  the  initial  endowments  satisfy  ea  G  int  and  the  u va 
continuously  converge  to  ua  then  Wu  lop-converge  ancillary  tightly  to  W. 

8.4  Theorem  Suppose  the  economies  { 8V ,  v  G  IV }  converge  to  8  and  that  the  Walrasians  Wu  as¬ 
sociated  with  them  are  Ky  Fan  bifunctions,  then  for  every  8V  comes  with  Walras  equilibrium  prices 
{pu  ,u  G  IV}.  Moreover,  any  sequence  of  such  equilibrium  points  has  at  least  one  cluster  point  and  any 
such  cluster  point  is  a  Walras  equilibrium  for  8. 

Proof.  As  the  sequence  Su  converges  to  8  and  the  functions  Wv  are  finite  on  A  x  A,  Theorem  8.1 
allows  us  to  infer  the  existence  of  equilibrium  points  for  each  one  of  the  economies  8U .  The  conclusions 
follow  from  Theorem  3.2  about  the  convergence  of  maxinf-points  and  the  fact  that  A  is  compact.  □ 


9  Convergence  of  solutions  to  generalized  equations 

There  is  an  extensive  literature  dealing  with  the  local  behavior  of  the  solutions  to  generalized  equations, 
or  inclusions,  under  perturbations  of  some  specific  parameters,  cf.  the  recent  monograph  of  Dontchev 
and  Rockafellar  [16]  and  references  therein;  see  also  the  work  of  Ait  Mansour,  in  particular  [1,  2].  A 
more  global  approach  was  already  followed  Aubin  and  Wets  [8]  but  with  a  boundedness  condition  on 
the  coderivatives  of  the  mappings.  Here,  the  conditions  imposed  are  on  the  properties  of  the  mappings 
themselves  rather  than  on  their  coderivatives.  Moreover,  we  limit  ourselves  to  a  simplified  situation; 
a  full  analysis  deserve  an  independent  treatment,  a  special  case  was  dealt  with  in  §5  as  an  alternative 
approach  to  the  convergence  of  variational  inequalities,  see  the  examples  that  follows  Proposition  5.1. 
Our  immediate  aim  is  to  illustrate  how  lop-convergence  can  be  exploited  to  get  us  on  the  way  to  con¬ 
vergence  results  for  the  solutions  of  Sv{x )  3  dv  to  the  solutions  of  S(x)  3  0  when  the  mappings  Su 
‘approximate’  S  and  dv  — »  0. 

Let  Su,  S  :  lRn  ^  Mm  be  convex- valued,  osc  mappings  and  consider  the  inclusions  Sv{x )  3  dv ,  S{x)  3  0 
with  du  -3  0.  For  x  G  domS'*'  =  {x  |  Su(x)  /  0},  let’s  designate  by  au(x ,  •)  the  support  function  of 
Su(x)  and  define  a(x,  •)  as  the  support  function  of  S.  The  bifunctions  that  get  associated  with  the  limit 
and  the  approximating  problems  are  defined  by 

K(x,  v)  =  a(x,  v)  and  Ku(x,  v)  =  (Tv(x,  v)  —  (du,  v), 

Again,  we  plan  to  rely  on  Ky  Fan’s  Inequality  for  existence  of  maxinf-points  of  these  bifunctions.  Of 
course,  for  all  x,  these  functions  are  convex  in  u.  The  need  for  the  upper  semicontinuity  of  these  support 
functions  leads  us  to  impose  some  further  conditions  on  the  collection  of  mappings  {S,SU,  u  G  IV }. 
Although  it’s  far  from  a  vital  condition,  let’s  carry  on  with  the  assumption  that  our  mapping  are 
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also  locally  bounded.  When  that’s  the  case,  the  outer  semicontinuity  of  the  mappings  is  enough  to 
guarantee  that  the  mappings  au(-,v),  cr(-,v)  are  use.  The  argument  is  provided  for  the  function  a,  it’s 
similar  for  av .  Indeed,  first  note  that  then  the  mapping  S  is  compact-valued  and  for  all  x  £  domS, 
sup  {(u,ii)  |  u  £  S'(x)}  is  attained  at  some  u  €  Su(x).  If  xk  — >  x  in  domS1,  the  corresponding  points 
uk  £  argmax  { (v,  u)  |  u  £  5(xfc)}  have  a  cluster  point  u  £  S(x)  as  follows  from  local  boundedness  and 
outer  semicontinuity  of  S.  Hence, 

limsupfc  a(xk,  v)  =  limsup k(v,uk)  <  o(x,v). 

Let’s  now  proceed  to  show  that  the  bifunctions  Ku  lop-converge  to  K  =  a  when  the  mappings  Sv 
graphically  converge  to  S,  i.e.,  gphS1'  — >■  gph5. 


9.1  Lemma  (range  of  a  compact- valued  osc  mapping).  Compact-valued  osc  mappings  defined  on  a 
compact  domain  have  compact  range  and,  hence,  are  locally  bounded. 

Proof.  Let  S  :  D  ^4-  Mm  be  such  a  mapping  and  let’s  show  that  rgei?  is  bounded.  To  the  contrary, 
suppose  rgeS  was  unbounded  which  means  there  exists  uk  £  S(xk)  with  \uk\/'oo  and  xk  £  D.  Since 
D  is  compact,  passing  to  a  subsequence,  if  necessary,  xk  — >  x  £  D  and  since  S  is  compact-valued  and 
osc  at  x  it  follows  that  for  any  e  >  0  arbitrarily  small  and  k  sufficiently  large,  S(xk)  C  S(£,x)  +  elB 
which  negates  the  possibility  of  having  such  a  sequence  {ufc}  and  consequently,  rgeS  must  be  bounded. 
That  rgeS  is  also  closed  follows  from  a  similar  argument:  when  uk  £  S(xk )  — »•  u,  again  passing  to  a 
subsequence  if  necessary,  implies  that  xk  — >  x  £  D  and  since  S  is  osc,  u  £  S(x)  C  rg eS.  □ 

9.2  Proposition  (inclusions:  lop-convergence  of  support  functions).  Suppose  du  —>■  0,  the  osc  (outer 
semicontinuous)  mappings  {S'1'  :  Dv  C  lRn  ^  Mm,  v  £  IV}  are  uniformly  locally  bounded  and  converge 
graphically  to  S  :  D  C  Mn  ^  dRm  ■  Then,  the  bifunctions 

{Ku  =  ou  -  (du,  •)  :  Du  x  Mm  -f  1R}  lop-converge  to  K  =  a  :  D  x  Mm  ->  ]R. 

Proof.  It  really  suffices  to  prove  that  the  support  functions  av  lop-converge  to  cr.  The  graphical 
convergence  of  the  mappings  Su  to  S  implies  that  gph5  is  closed  and  thus,  S  is  also  osc.  Moreover, 
since  the  graphical  limit  at  some  point  x  £  D  can  be  expressed  as 

[J  Limsup, Su(xv)  C  S(x)  C  [J  Liminf^oo  Su(xu), 

{xv— >x}  {xv^x} 

see  [33,  Proposition  5.33],  it  immediately  follows  that  S  is  locally  bounded  since  the  mappings  Su  are 
uniformly  locally  bounded. 

Let’s  now  turn  to  condition  3.1(a).  Given  v  and  (xu  £  Du )  — >•  x,  we  have  to  exhibit  vu  — >  v 
such  that  limsup^  ou(xv,  vu)  <  a(x,v).  Let’s  simply  choose  vu  =  v  and  let  uv  £  argmax  { (v,  it)}  |  u  £ 
Su(x v)}.  Now  recall  that  the  mappings  Su  and  S  are  uniformly  locally  bounded  and,  consequently, 
every  subsequence  of  {uu ,  v  £  LV}  comes  with  a  further  converging  subsequence.  For  any  cluster  point 
u  of  (uu ,  v  £  IN},  restricting  our  attention  to  the  subsequence  converging  to  uu  fi,  N  C  IN,  one  has 

lirn u£]y  ou (xu ,v)  =  lim^eAr(uy,  v)  =  (u,v)  <  a(x,v) 
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from  which  one  immediately  concludes  that  limsup  ou(xu,  v)  <  a(x,v). 

Condition  3.1(b)  is  verified  as  follows.  For  x  £  D  and  vv  — >  v,  let  u  £  argrnax  { (v,  u)  |  u  £  S'(x)}. 
Graphical  convergence,  in  particular  Liminf,,  gph  Su  D  gph  S,  implies  that  there  exists  a  sequence 
(xv ,uu)  £  gph  S'1'  — >  (x,u).  Taking  this  into  account,  with 

(i)  o(xu,vu)  >  (vu,uu),  and 

(ii)  (vv,uv)  ->  (v,u), 

one  obtains,  liminf uo(xu,vu)  >  a(x,v)  as  required.  □ 


10  Set-valued  fixed  points 

Given  C  a  compact,  convex  subset  of  lRn  and  S  :  C  ^  C  C  lRn,  a  convex-valued  osc  mapping, 
Kakutani’s  extension  of  Brouwer  fixed  point  problem  deals  with 

finding  x  £  C  such  that  x  £  S(x). 

It’s  easy  to  verify  that  a  solution  of  such  a  problem  is  also  a  maxinf-point  of  the  mapping  F  :  CxC  — »  1R 
with 

F(x,  v )  =  sup  { (x  —  v,  z  —  x)  |  z  £  S(x)  Cd} 
and  vice-versa.  Our  approximating  problems  and  the  associated  bifunctions: 

find  x  £  Cv  such  that  x  £  Su(x), 

Fu(x,  v )  =  sup  {(x  —  v,  z  —  x)  |  2  €  Sv{x)  C  C1'} 

with  the  finite- valued  bifunctions  Fu  defined  on  Cv  x  Cv .  Existence  issues  can  again  be  resolved  via 
an  appeal  to  Ky  Fan’s  inequality  since  the  bifunctions  F,  Fu  are  clearly  convex  in  v,  non-negative  on 
the  diagonal,  i.e. ,  F(x,x)  >  0  on  C  x  C  and  Fu(x,x)  >  0  on  Cv  x  Cv ,  and  upper  semicontinuous  in 
x  as  follows  immediately  from  the  properties  of  sup-projection  [38]  when  the  mappings  S,  Sv  are  outer 
semicontinuous.  Convergence  of  the  maxinf-points  and  the  approximate  solutions  (e-maxinf-points) 
will  follow  from  the  results  in  §3,  if  we  come  up  with  conditions  under  which  the  lop-convergence  of  the 
bifunctions  Fu  to  F  is,  not  just  ancillary  tight,  but  actually  tight. 

10.1  Theorem  (convergence  of  the  fixed  points  of  set-valued  mappings).  Consider  the  (set-valued) 
fixed  points  problems  described  above  and  their  associated  bifunctions  F  and  Fu .  If  the  sets  Cu  — >  C 
and  the  osc  mappings  Sv  converge  graphically  to  S,  i.e.,  gph  Sv  gph  S,  then  the  bifunctions  Fu 
lop-converge  tightly  to  F.  Consequently,  for  any  £u\0,  any  cluster  point  of  a  sequence  of  ev-maxinf 
points  is  a  maxinf-point  of  F  and  every  fixed  point  of  the  limit  problem  is  the  limit  of  approximate 
fixed  points  of  the  approximating  problems. 

Proof.  Since  C  is  compact  and  convex,  it  follows  that  for  v  sufficiently  large,  the  sets  Cu  must  lie  in 
a  bounded  region  [33,  Corollary  4.12].  We  proceed  as  if  this  is  the  case  for  all  u.  This  implies,  Lemma 
9.1,  that  the  ranges  of  the  mappings  S  and  Su  are  uniformly  bounded.  Hence,  verifying  3.2(a-t)  and 
3.2(b-t)  really  comes  down  to  checking  3.2(a)  and  3.2(b). 
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For  any  sequences  vv  £  Cv  — >•  v  £  C  and  xv  £  Cv  ^  x  £  C,  since  gph  Su  — >  gph  S,  the  functions 
gv(z)  =  (xu  —  vu,  z  —  xv)  defined  on  Sv{xv)  C  Cv 

hypo-converge  tightly  to 

g(z)  =  (x  —  v,  z  —  x)  defined  on  S(x)  C  C. 

We  can  apply  [25,  Theorem  2]  and  with  zv  £  argma xgu  — recall  that  Sv(xu)  is  compact — ,  one  has 
gv(zu)  — >•  g{z)  =  F(x,v)  where  z  is  any  cluster  point  of  the  sequence  {zu}u&n.  Whence,  both  3.2(a) 
and  3.2(b)  are  satisfied  now  implying  that  the  bifunctions  Fu  lop-converge  tightly  to  F  and  the  full 
implications  of  Theorem  3.2,  Theorem  3.3  and  its  Corollary  3.4  are  available.  □ 

10.2  Example  (fattening  up  a  function  with  a  fixed  point).  This  problem  was  brought  to  our  attention 
by  Jong-Shi  Pang  in  connection  with  his  work,  and  that  of  his  collaborators,  on  cognitive  radio  games 
[29,  34,  18].  Let  F  :  C  -»  C  be  continuous  on  C  a  nonempty  convex  compact  subset  of  IR'1  such  a 
mapping  has,  at  least,  one  fixed  point,  say  x  £  C,  Theorem  6.1.  What  are  the  the  conditions  under 
which  one  can  ‘ fatten  up ’  F,  say  F£  :  C  ,  such  that  their  fixed  points  as  e\0  will  converge  to  those 
of  F,  and  in  particular,  to  x. 

Detail.  It’s  easy  to  see  that  one  can’t  proceed  arbitrarily.  A  minimal  requirement  is  to  make  sure 
that  these  mappings  { F£  :  C  lRn,  e\0}  are  outer  semicontinuous  (osc)  in  order  to  guarantee  the 
existence  of  fixed  points,  as  well  as  approximate  fixed  points.  One  might  be  tempted  to  require,  in  fact, 
that  these  mappings  are  continuous,  i.e. ,  are  also  inner  semicontinuous,  so  as  to  be  able  to  appeal  to  a 
‘generalized’  implicit  function  theorem.  However,  this  a  pretty  strong  requirement,  not  well  suited  to 
the  applications  [29].  Let’s  take  for  granted  that  are  we  only  referring  to  0  <  e  <  e  <  oo.  Then,  these 
mappings  {.Fe}£>0  are  uniformly  locally  bounded,  in  fact,  their  range  is  (uniformly)  bounded.  This 
brings  us  in  the  reach  of  the  previous  theorem  which  tells  us  that  if  the  fattened  up  mapping  graphically 
converge  to  F  as  e\0,  we  are  assured  that  any  sequence  of  fixed  points  of  the  mappings  F£  will  cluster 
to  a  fixed  point  of  F,  Moreover,  via  Corollary  3.4,  one  can  always  find  a  sequence  of  approximate  fixed 
points  of  the  mappings  F£  that  converges  to  a:.  □ 


References 

[1]  M.  Ait  Mansour.  Sensitivite  et  stability  des  points  d’equilibre  et  quasi- equilibre:  applications  aux 
inegalites  hemi-variationelles  et  variationelles.  PhD  thesis,  Universite  Cadi  Ayyad,  Marrakech, 
2002. 

[2]  M.  Ait  Mansour  and  H.  Riahi.  Sensitivity  analysis  for  abstract  equilibrium  problems.  Journal  of 
Mathematical  Analysis  and  Applications,  306:684-691,  2005. 

[3]  K.  Arrow  and  G.  Debreu.  Existence  of  an  equilibrium  for  a  competitive  economy.  Econometrica, 
22:265-290,  1954. 

[4]  H.  Attouch,  D.  Aze,  and  R.J-B  Wets.  Convergence  of  convex-concave  saddle  functions:  continu¬ 
ity  properties  of  the  Legendre-Fenchel  transform  with  applications  to  convex  programming  and 
mechanics.  Annales  de  I’Institut  H.  Poincare:  Analyse  Nonlineaire,  5:537-572,  1988. 


27 


[5]  H.  Attouch  and  R.J-B  Wets.  Convergence  des  points  min/sup  et  de  points  fixes.  Comptes  Rendus 
de  I’Academie  des  Sciences  de  Paris ,  296:657-660,  1983. 

[6]  H.  Attouch  and  R.J-B  Wets.  A  convergence  theory  for  saddle  functions.  Transactions  of  the 
American  Mathematical  Society,  280:1-41,  1983. 

[7]  J.-P.  Aubin  and  I.  Ekeland.  Applied  Nonlinear  Analysis.  J.  Wiley  Intersciences,  1984. 

[8]  J.-P.  Aubin  and  R.J-B  Wets.  Stable  approximations  of  set-valued  maps.  Annales  de  I’Institut  H. 
Poincare:,  serie  C,  5:519-535,  1988. 

[9]  A.  Bagh.  Epi/hypo-convergence:  the  slice  topology  and  saddle  points  approximation.  Journal  of 
Applied  Analysis,  4:13-39,  1996. 

[10]  X.  Chen,  R.J-B  Wets,  and  Y.  Zhang.  Stochastic  variational  inequalities:  Residual  minimization, 
smoothing/sample  average  approximation.  Technical  report,  The  Hong  Kong  Polytechnic  Univer¬ 
sity,  2011. 

[11]  R.W.  Cottle,  J.-S.  Pang,  and  R.E.  Stone.  The  Linear  Complementarity  Problem.  Academic  Press, 
1992. 

[12]  J.  Danskin.  The  Theory  of  Max-Min.  Springer,  1967. 

[13]  G.  Debreu.  Theory  of  Value.  J.  Wiley  Sz  Sons,  1959. 

[14]  J.  Deride.  Computing  and  sensitivity  analysis  for  an  computing  and  sensitivity  analysis  for  an 
equilibrium  in  incomplete  markets:  A  duality  approach.  Master’s  thesis,  Universidad  de  Chile, 
2011. 

[15]  J.  Deride,  A.  Jofre,  and  R.J-B  Wets.  Computing  equiibrium  points  in  incomplete  markets:  a 
variational  approach.  In  2010-ICCOPT  Meeting,  Santiago,  2010. 

[16]  A.  Dontchev  and  R.T.  Rockafellar.  Implicit  Functions  and  Solution  Mappings.  Springer  Mono¬ 
graphs  in  Mathematics.  Springer,  2009. 

[17]  F.  Facchinei  and  J.-S.  Pang.  Finite- Dimensional  Variational  Inequalities  and  Complementarity 
Problems.  Springer,  2003. 

[18]  F.  Facchinei  and  J.-S.  Pang.  Convex  Optimization  in  Signal  Processing  and  Communications. 
Cambridge  University  Press,  2009. 

[19]  K.  Fan.  A  minimax  inequality  and  applications.  In  O.  Shisha,  editor,  Inequalities  -  III,  pages 
103-113.  Academic  Press,  1972. 

[20]  F.  Flores-Barzan  and  R.  Lopez.  The  linear  complementarity  problem  under  asymptotic  analysis. 
Mathematics  of  Operations  Research,  30:73-90,  2006. 

[21]  S.  Gowda  and  J.S.  Pang.  The  basic  theorem  of  complementarity  revisited.  Mathematical  Program¬ 
ming,  58:161-177,  1993. 


28 


[22]  G.  Giirkan,  A.  Ozge,  and  S.  Robinson.  Sample-path  solution  of  stochastic  variational  inequalities. 
Mathematical  Programming,  84:313-333,  1999. 

[23]  G.  Giirkan  and  J.-S.  Pang.  Approximations  of  Nash  equilibria.  Mathematical  Programming,  Ser. 
B,  117:223-253,  2009. 

[24]  P.T.  Harker  and  J.-S.  Pang.  Finite  dimensional  variational  inequality  and  nonlinear  complemen¬ 
tarity  problems:  A  survey  of  theory,  algorithms  and  applications.  Mathematical  Programming, 
48:161-220,  1990. 

[25]  A.  Jofre  and  R.  Wets.  Variational  convergence  of  bivariate  functions:  lopsided  convergence.  Math¬ 
ematical  Programming,  Ser.  B,  116:275-295,  2009. 

[26]  A.  Jofre  and  R.J-B  Wets.  Continuity  properties  of  Walras  equilibrium  points.  Annals  of  Operations 
Research,  114:229-243,  2002. 

[27]  V.I.  Konnov.  Equilibirum  Models  and  Variational  Inequalities,  volume  210  of  Mathematical  Science 
and  Engineering.  Elsevier,  2007. 

[28]  D.T.  Luc  and  R.J-B  Wets.  Outer  semicontinuity  of  positive  hull  mappings  with  application  to 
semi-infinite  and  stochastic  programming.  SIAM  Journal  on  Optimization,  19:700-713,  2008. 

[29]  J.-S.  Pang  and  G.  Scutari.  Nonconvex  cognitive  radio  games-part  i:  quasi-Nash  equilibria,  (sub¬ 
mitted  to  IEEE  Transactions  on  Signal  Processing),  University  of  Illinois,  Industrial  and  Enterprise 
Systems  Engineering,  2010. 

[30]  E.  Polak.  Optimization  Algorithms  and  Consistent  Approximations.  Springer,  1997. 

[31]  S.M.  Robinson.  Analysis  of  sample-path  optimization.  Mathematics  of  Operations  Research, 
21:513-528,  1996. 

[32]  R.T.  Rockafellar.  Convex  Analysis.  Princeton  University  Press,  1970. 

[33]  R.T.  Rockafellar  and  R.J-B  Wets.  Variational  Analysis.  Springer,  1998  (3rd  printing  2009). 

[34]  G.  Scutari  and  J.-S.  Pang.  Joint  sensing  and  power  allocation  in  nonconvex  cognitive  radio  games: 
Nash  equilibria  and  distributed  algorithms,  (submitted  to  IEEE  Transactions  on  Information 
Theory),  SUNY  Buffalo,  Electrical  Engineering,  2011. 

[35]  D.  Walkup  and  R.J-B  Wets.  Continuity  of  some  convex-cone  valued  mappings.  Proceeding  of  the 
American  Mathematical  Society,  18:229-235,  1967. 

[36]  L.  Walras.  Elements  d’economie  pure;  ou,  Theorie  de  la  richesse  sociale.  L.  Corbaz  &  cie,  1874. 

[37]  R.J-B  Wets.  On  the  continuity  of  the  value  of  a  linear  program  and  of  related  polyhedral-valued 
multifunctions.  Mathematical  Programming  Study,  24:14-29,  1985. 

[38]  R.J-B  Wets.  Lipschitz  continuity  of  inf-projections.  Computational  Optimization  and  Applications, 
25:269-282,  2003. 


29 


